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I. INTRODUCTION 

In 1973, Altschuler and Frantz patented an idea for 
creating an atom interferometer (Altschuler and Frantz, 
1973). The beam splitter in their apparatus was a 
standing-wave optical field. Their ideas were rekindled 
by Dubetsky et al. (1984), who presented detailed cal- 
culations of atomic scattering by standing-wave fields in 
the context of atom interferometry. It was not until re- 
cently, however, that experimentalists were successful in 
constructing the first atom interferometers. Double-slit 
interference (Carnal and Mlynek, 1991; Shimizu et al., 
1992), Fraunhofer diffraction by microfabricated struc- 
tures (MS) (Keith et al., 1991; Ekstrom et al., 1995) 
or by resonant standing wave fields (SW) (Rasel et al., 
1995; Giltner et al., 1995), and Fresnel diffraction by one 
(Chapman et al., 1995) or two (Clauser and Li, 1994) 
MS have all been observed using atomic beams as mat- 
ter waves. 

Two types of atom-optical elements have been used as 



beam splitters in these experiments, standing wave fields 
or microfabricated structures. A SW beam splitter al- 
lows one to operate with relatively dense atomic beams, 
having densities up to 10 10 cm -3 and flow densities up 
to 10 15 cm _2 s _1 . Moreover, by varying the atom field 
detuning, one can use SW beamsplitters as either am- 
plitude or phase gratings. Additional degrees of freedom 
are provided by the polarization of the field which can act 
selectively on targeted magnetic state sublevels. A the- 
ory of atom interference in standing wave fields has been 
developed by Altschuler and Frantz (1973) , Dubetsky 
et al. (1984), Chebotayev et al. (1985), Borde (1989), 
Friedberg and Hartmann (1993,1993a), Dubetsky and 
Berman (1994) and Janicke and Wilkens (1994). In con- 
trast to SW beam splitters, MS usually scatter atoms in 
a state independent manner; as a consequence, most ex- 
periments involving MS use atoms in their ground (or, 
possibly, metastable) states. Microfabricated structures 
provide 100% modulation of the incident atomic beam. 
They offer the additional advantage that their period and 
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duty cycle (ratio of slit opening to period) can be cho- 
sen arbitrarily within the limits of current lithographic 
technology. A theory of atom interference using MS has 
been developed by Turchette et al. (1992) , Clauser and 
Reinsch (1992) and Carnal et al. (1995) . 

Both the splitting of an atomic beam into two or more 
beams coherent with respect to one another and the re- 
combining of the scattered beams are physical processes 
that are essential to the operation of an atom interfer- 
ometer. We consider scattering of atoms by an ideal MS, 
having an infinite number of slits, period d, duty cycle 
/, and 100% transmission through the slits. Each MS is 
normal to the y-axis, and the slits are oriented in the z- 
direction, so that the axis of the MS is in the x-direction 
(see Fig. @). 
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FIG. 1. Matter- wave interference (a) and shadow effect (b) 
for a collimated atomic beam incident on a microfabricated 
structure (MS) . (a) When an atom having given center of 
mass momentum p || y scatters from a MS having wave vec- 
tor k =2nit/d, where d is the period of the MS, the output 
wave function consists of a set of momenta states p+nhk 
(two of them |p) and |p+?ik} are shown). A superposition 
of these states' wave functions leads to an interference pat- 
tern on the screen, having the same period d as the MS's. 
(b) When a collimated beam of particles moving along classi- 
cal trajectories scatters from MS, the beam profile imprinted 
by MS is copied on the screen. A beam consisting of two 
velocity subgoups, having zero (v = 0) and nonzero (v ^ 0) 
velocity projections on the a;— axis is shown. Gratings associ- 
ated with these two subgroups are shifted on the screen from 
on another by Ax = vt — JL, where t is the time of flight 
between the MS and the screen. 

After scattering from a MS, each in-coming atomic 
state i\) having x-component of center-of-mass momen- 
tum p splits into a set of out-going states ip n having x- 
components of momenta p + nhk (k — 2it/d, n is an 
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(1) 



Interference of two components, such as "0o aud ipi , on a 
screen (see Fig. 0a) leads to an atomic density grating 



p oc Re [-!/>() ^i] = cos [kx] 



(2) 



having the same period d as the MS. 

Observation of this grating in the experiments listed 
above often has been considered as direct evidence for 
matter- wave interference. Nevertheless, one can easily 
see that such a conclusion is not necessarily justified. For 
particles moving along classical trajectories (Fig. |l|b), 
and for a beam whose angular divergence is sufficiently 
small to satisfy 



(3) 



where L is a distance on the order of the distance be- 
tween the MS and the screen, a shadow of the MS can 
be seen on the screen at distances L where all matter- 
wave effects are completely negligible. This example is 
the simplest manifestation of the classical shadow effect 
(Chebotayev et al., 1985; Dubetsky and Berman 1994). 

To distinguish quantum matter-wave interference from 
the classical shadow effect, one needs to observe addi- 
tional features of the phenomena or to choose a scheme 
where one of the effects is excluded. Young's double-slit 
experiment (Carnal and Mlynek, 1991; Shimizu et al., 
1992), as well as interference produced by a phase grat- 
ing created using light that is far-detuned from atomic 
transition frequencies (Rasel et al., 1995), cannot be ex- 
plained in terms of atoms moving on classical trajecto- 
ries. A matter-wave interpretation is also necessary if 
one observes a shift in the fringe pattern resulting from 
an index change in one of the arms of an interferometer 
(Ekstrom et al., 1995). We determine below those partic- 
ular conditions for which pure quantum interference can 
be obtained using MS. 

Let us estimate a typical distance for which quantum 
interference effects have to be included. Consider an in- 
cident beam which has zero angular divergence. The 
atoms are assumed to be in a pure state having mo- 
mentum p = (0,p y ,0) (see Fig. 0). Localization of the 
atoms inside each slit leads to an uncertainty in the x- 
component of atomic momentum 5p ~ § , where it is as- 
sumed that the slit width fd is comparable with the MS 
period d. A beam passing through the slits acquires an 
angular divergence 59 ~ — , and atoms passing through 

Py 

a given slit are deposited on the screen with a spot size 
Sx ~ L89 ~ L-fi-, Interference occurs when spots pro- 

ap v 

duced by neighboring slits overlap, i. e. when Sx ~ d. 
One finds, therefore, that the characteristic distance for 



2 



which matter-wave interference plays an essential role is 
given by 



L~ Ln 



(4) 



where 



L T = 2d 2 /\ dB (5) 



is the so-called Talbot distance and XdB = — is the 

Py 



atomic dc Broglie wavelength. 
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FIG. 2. The incident atomic wave function after scatter- 
ing from MS transforms into a superposition of the divergent 
waves emitted from the each slit, having angular divergence 
89 inverse to the MS's period d. (for the given duty cycle). 
One estimates Talbot distance as a distance between MS and 
screen large enough to provide overlapping of neighboring di- 
vergent waves and their interference. 

The manner in which this distance appears in the the- 
ory of the optical or atomic Talbot-effect is well known 
in the context of the Fresnel-Kirchhoff theory of diffrac- 
tion [see, for example, (Patorski, 1989; Winthrop and 
Worthington, 1965) for optical Talbot-effect theory and 
(Chapman et ah, 1995; Turchette et al., 1992; Clauser 
and Reinsch, 1992; Carnal et al., 1995) for atomic Talbot- 
effect theory]. One can obtain Eq. (Q, ||) using another 
description. When the number of slits in the MS and the 
area of the incident atomic beam are infinite, construc- 
tive interference occurs only for those directions in which 
p is changed by an integral multiple of the recoil mo- 
mentum hk. In the atomic "rest- frame" (a frame moving 
along the y— axis with velocity 



: p y /M, 



(6) 



where M is the atomic mass) an out-going state with mo- 
mentum p + nhk ([!]) acquires a phase <j) n = e p + n rikt/fo, 

2 

where e p = is the kinetic energy associated with 
atomic motion along the ir-axis and t is the time after 
scattering from the MS. Comparing this phase with the 
phase <po that an atom would acquire in the absence of 
the MS, one sees that a dephasing 5<f) = <fr n — 4>o occurs 



for the different out-going state amplitudes as a result of 
diffraction. The relative dephasing is 



where 



n4> D + n <j) t , 



bp = kvt, 



4>t = Wfct, 

v = -ft* is the x component of atomic velocity, and 



hk 2 

2M 



(7) 

(7a) 
(7b) 

(8) 



is a recoil frequency, related to the energy that an atom, 
having initial momentum p = 0, acquires as a result of 
recoil during scattering. The two contributions to the 
dephasing (ffi) have different origins. The phase <f>o is a 
Doppler phase that does not disappear in the classical 
limit h — > 0; consequently, it must be classical in nature. 
It is evident from Fig. [l]b that, for atoms incident with a 
nonzero value of p, the shadow moves along the x— axis 
together with atoms. The nodes in the shadow are dis- 
placed by a distance Ax — vt along the a;— axis, which 
corresponds to a phase shift of 27r-^ for the atomic grat- 
ing, a phase shift that coincides with <f>jj . The phase 4>r> 
is analogous to the phase a moving dipole driven by an 
optical field would acquire in its rest frame as a result of 
the Doppler frequency shift. 

The phase 4>t hi Eq- (0) is a quantum addition to 
the dephasing, resulting from recoil. This contribution 
is responsible for atomic scattering and for matter-wave 
interference. Quantum effects have to be included when 
this phase is of order unity, that is, for times of order 



t T = 2-K/lxlk- 



(9) 



One finds that length associated with this time in labo- 
ratory frame coincides with the Talbot distance 



Lt = uIt- 



(10) 



Thus, we are led to the same conclusion that we reached 
above by considering the scattering from two adjacent 
slits; for cuktr 1 or L > Lt, one must use a quantized 
description of the atomic center-of-mass motion. 

In this chapter we consider the Talbot-effect and other 
interference phenomena as a consequence of the recoil- 
effect. In the context of the nonlinear interaction of opti- 
cal fields with an atomic vapor, the recoil effect was con- 
sidered by Kol'chenko et al. (1968) and observed by Hall 
et al. (1976). Quantum structure resulting from the scat- 
tering of atoms by a resonant standing wave (resonant 
Kapitza-Dirac effect), which can be attributed to atomic 
recoil, was discussed theoretically by Kazantsev et al. 
(1980) and observed by Moskovitz et al. (1983). Splitting 
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of optical Ramsey fringes (Baklanov et al., 1976) associ- 
ated with the resonant Kapitza-Dirac effect was discussed 
theoretically by Dubetsky and Semibalamut (1978) and 
observed by Barger et al. (1979). Matter-wave interfer- 
ence resulting from resonant Kapitsa-Dirac scattering in 
a standing wave field has been studied theoretically by 
Altschuler and Frantz (1973), Dubetsky et al. (1984) 
and observed by Rasel et al. (1995). The theory of atom 
interference presented here, based on an interpretation of 
scattering of atoms by MS in terms of the recoil effect, is 
a natural extension of the work involving standing-wave 
fields. 

This chapter is organized as follows: in the next Sec- 
tion we discuss conditions necessary for the observation 
of matter-wave interference in different regimes. Rigor- 
ous proof of the equivalence of theories based on Fresnel- 
Kirchhoff integrals and on the recoil effect is given in Sec- 
tion III, as is a discusion of the atomic gratings that can 
be produced as a consequence of the Talbot effect. The 
classical shadow effect is analyzed in Section IV. Section 
V is devoted to a theory of the Talbot-Lau effect. The 
Talbot and Talbot-Lau effects for a thermal beam are 
considered in the Section VI, in the limit where the char- 
acteristic length scale in the problem is larger than the 
Talbot length Lt- A discussion of the results is given in 
Section VII. 



II. QUALITATIVE CONSIDERATIONS 

The scattering of atoms by gratings can be separated 
roughly into three categories: classical scattering, Fres- 
nel diffraction, and Fraunhofer diffraction. The limit of 
Bragg scattering (Martin et al., 1988), in which uiki' /u > 
1, where £' is the grating thickness, is not discussed in 
this chapter. 



A. Fresnel Diffraction 

The Fresnel diffraction limit occurs when Ukt ~ 1 or 
L ~ Lt- Owing to the angular divergence of the inci- 
dent beam, it is possible that the diffraction pattern at 
L ~ Lt will be washed out. To ensure that this does not 
occur, it is necessary that the spread of Doppler phases, 
kutTOb — kLT0b, be smaller than unity. This require- 
ment corresponds to inequality (0) when L ~ Lt- Using 
Eq. (pf), the condition on 9b can be restated as 



where 



Qb < 



Tik 
Py 



(11) 



(12) 



is the angle associated with a single atomic recoil at the 
MS (Sp = hk). The Talbot effect refers to the self- 
imaging of a grating in the Fresnel diffraction limit. For 
self-imaging to occur, the displacement of the scattered 
atomic beam LtOj, must be much smaller than the beam 
diameter D, which translates into the condition 



(13) 



In this limit, one can consider the beam diameter to be 
infinite to first approximation; finite beam effects (or, 
equivalently, gratings with finite slit number) are dis- 
cussed by Clauser and Reinsch (1992) . 

Conditions (g, [□} [jj) are sufficient to observe the Tal- 
bot effect. In this case, the contribution to the wave 
function's phase resulting from atomic recoil is of order 
unity, the scattered beams overlap almost entirely with 
one another on the screen, and the atomic gratings are 
not washed out after averaging over atomic velocities v in 
the incident beam. Matter-wave interference is a critical 
component of the Talbot effect. 



B. Fraunhofer Diffraction 

Although the Talbot effect illustrates a matter-wave 
interference phenomenon, it does not result in an atom 
interferometer having two arms that are nonovcrlapping. 
We refer to the Fraunhofer diffraction limit as one in 
which the various diffraction orders are nonoverlapping 
at a distance L from a single grating. The grating then 
serves as a beam splitter that physically separates the 
incident beam into two or more beams. To physically 
separate the various diffraction orders over a distance L, 
one must require that 



L8 d > D. 



(14) 



Using the fact that 9d = XdB/d and setting t = ^, one 
can recast this inequality as 



uj k t > D/d 



L/Lt > D/d, 



(15) 



which requires the quantum phase 4>t (7b) to be larger 
than D/d >> 1. Consequently, quantum effects play an 
essential role in an atom interferometer having nonover- 
lapping beam paths. Note that the Fraunhofer limit can- 
not be reached for a beam having infinite diameter. We 
do not consider matter-wave interference in the Fraun- 
hofer limit in this work. 



C. Talbot-Lau Regime. 

To achieve spatial separation of the beams in the 
Fraunhofer limit and to observe the Talbot effect, the 
angular divergence 9b of the incident beam must be less 
than 9d- For typical values d ~ 200 nm, M ~ 20 



4 



A.u., u ~ 10 5 cm/s, the deflection angle Od ~ 10~ 4 
rac?. Atomic beams having 9 b <C 6*d have been used to 
observe atomic scattering by standing waves (Moskovitz 
et. al., 1983), to build a two-arm atom interferometer 
(Keith et al., 1991; Rasel et al, 1995; Giltner et al, 
1995), and to observe the Talbot-effect (Chapman et al., 
1995). Such strong collimation results in a decrease in 
the atomic flux and a corresponding decrease in signal 
strength that may be a limiting factor in certain applica- 
tions of matter-wave interference, such as atom lithogra- 
phy (Timp et al., 1992). Alternatively, one can observe 
matter-wave interference in beams having larger angular 
divergence using the atomic Talbot-Lau effect [see for ex- 
ample (Patorski, 1989)]. In the atomic Talbot-Lau effect 
two or more MS are used. Doppler dephasing following 
the first MS washes out the normal Talbot effect, but 
subsequent scattering by a second MS can result in a 
Doppler rephasing that ultimately leads to a Talbot-likc 
interference pattern. The dephasing-rephasing process is 
analogous to that occurring in the production of pho- 
ton echoes (Dubetsky et al., 1984). The atomic Talbot- 
Lau effect has been observed recently by Clauser and Li 
(1994) using a K beam scattered by MS. The incident 
beam is not separated into nonoverlapping beams in the 
Talbot-Lau effect, but the origin of the interference pat- 
tern can still be traced to matter-wave interference since 
it is related to Fresnel diffraction. 

D. Classical Scattering 

It is worthwhile at this point to return to the classical 
shadow effect. The shadow effect in a collimated beam is 
obvious; the atomic grating produced by the MS simply 
propagates in space over a distance in which diffraction 
can be ignored (L/Lt C 1). If one has an ensemble of 
atoms incident on a MS from different angles, each veloc- 
ity subgroup creates its own grating. Just after passing 
through the MS, all the gratings are the same, having the 
profile of the MS. Downstream from the MS, the atomic 
gratings corresponding to different velocity classes move 
in different directions (two of them are shown in Fig.[l]b) . 
The grating of the ensemble as a whole is washed out at 
a distance £ from the MS, provided that 

£6 b > d. (16) 

For 

L6 b > d, (17) 

the shadow of the atomic grating is washed out well be- 
fore reaching the screen or any subsequent MS in the 
experimental setup. 

There are two ways to view the washing out of the 
grating. As discussed above, the washing out is a result 
of the different classical trajectories of the atoms. An 



alternative view allows one to relate this phenomenon to 
that encountered in optical coherent transients. Imagine 
that the incident beam consists of a number of velocity 
subgroups, having v x = v ~ u6 b . When the atoms are 
scattered by a MS having a transmission function Xi( x )j 
the atomic density immediately following the grating is 
proportional to xi( x )- Downstream from the MS, the 
atomic grating in the ^-direction simply propagates with 
velocity v, leading to a density distribution that varies 
as Xi( x ~ v t)> where t — y/u and y is the distance from 
the grating. If this density distribution is expanded in a 
Fourier series, one finds terms in the sum that vary as 
cos[nk(x — vt)], where n is an integer. In this picture, 
particles in velocity subgroup v acquire a Doppler phase 
(f7a| ) of order 

4>d ~ kvt ~ O b L/d » 1, (18) 

as they propagate a distance of order L from the MS. 
The decay of the macroscopic grating is analogous to the 
free-induction-decay of the macroscopic polarization of a 
Doppler-broadened atomic vapor following excitation by 
an optical pulse. In this approach, one can draw on many 
processes that are well known in the theory of optical co- 
herent transients. 

Although the atomic grating is washed out following 
the interaction with the MS, it is possible to restore the 
original macroscopic atomic grating by placing a second 
MS between the first MS and the screen. This effect has 
been observed recently by Batelaan et al. (1996) using a 
metastable Ar beam. In the classical trajectory picture, 
the restoration of the grating corresponds to a Moire pat- 
tern. In the Doppler dephasing picture, the restoration 
is analogous to the dephasing-rephasing process that oc- 
curs for a photon echo. The first MS starts a dephasing 
process for the different velocity subgroups and the sec- 
ond MS results in a rephasing process (see below, section 
IV. A). At a particular focal plane, where the rephasing 
is complete, a macroscopic grating appears. For 

L <C Lt (19) 

effects relating to quantization of the atomic center-of- 
mass motion play no role jjj . 

If L ~ Lt, one has to include recoil effects. Gratings 
appearing in this regime are usually associated with the 
Talbot-Lau effect, i. e. with the interference of light for 
the optical case or quantum interference of matter waves 
in the case of the atomic Talbot-Lau effect. It is shown 
in Sec. V, however, that the position of the focal planes 
and gratings' periods are often the same as in the classi- 
cal case. From this point of view, the Talbot-Lau effect 
is a quantum generalization of the shadow effect. 

In summary, one can conclude that interference is 
qualitatively different for collimated beams and beams 
with large angular divergence. For collimated beams 
{9 b L <C d) one has two interesting regimes, 
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L~L T [u k t~ 1] (20) 

and 

L/L T > D/d [u k t > D/d] , (20a) 

corresponding to Fresnel (Talbot effect) and Fraunhofer 
diffraction (nonoverlapping scattered beams), respec- 
tively. For beams having angular divergence 9bL ^> d, 
the Fresnel and Fraunhofer diffraction patterns would 
wash out following a single MS. Restoration of the atomic 
gratings in this limit can be achieved using two or more 
MS. Distances L ~ Lt correspond to the atomic Talbot- 
Lau effect and distances 

L <C Lt w fe i<l (21) 

correspond to the classical shadow effect. 

III. TALBOT EFFECT 
A. Talbot effect as a recoil effect. 

In this section we show that the Talbot effect is a con- 
sequence of the recoil an atom undergoes when it passes 
through a microfabricated structure (MS). We assume 
that the MS is located in the plane y = 0, normal to 
the direction of propagation of the atomic beam. The 
MS consists of an infinite number of slits oriented in 
the z— direction; as such, only the a;— dependence of the 
atomic wave function changes when atoms pass through 
the slits. Atomic motion in the y— direction can be con- 
sidered as classical in nature provided that A<js /d <C 1, 
but motion along the a;-axis must be quantized. In this 
section, it is assumed that the incident beam is strongly 
collimated, Ob <C d/L, where d is the period of the MS 
and L is the distance from the MS to the screen. As such, 
we can neglect any spread in the transverse velocities in 
the initial beam and consider all atoms to be incident 
with transverse momentum p = 0. 

After passing through the MS, the wave function for 
an atom is given by 

i>{x) = r)(x), (22) 

where rj(x) is the amplitude transmission function asso- 
ciated with the MS. In the momentum representation, 
ij){x) can be written as a superposition of states having 
momenta p — mhk, where m is an integer and k = 2tt j d. 
Explicitly, one finds that the Fourier transform of ip(x) 
is given by 

4>(p) — V2irh Tj m 8{jp — mhk), (23) 

m 

where 



/dr 
—e-* mk * V (x) (24) 

is a Fourier coefficient. Unless indicated otherwise, all 
sums run from — oo to +oo. The terms with m =/= in 
Eq. ( p3| ) can be associated with atomic scattering at an- 
gles mhk/py, where p y is the longitudinal momentum in 
the atomic beam. It is assumed that p y is constant for all 
atoms in the beam - this restriction is relaxed in Sec. VI. 
In analogy with electron scattering from a standing wave 
field [Kapitza-Dirac effect (Kapitza and Dirac, 1933)] or 
atomic scattering from a resonant standing wave field 
[resonant Kapitza-Dirac effect (Kazantsev et al., 1980)], 
the scattering from the MS can be interpreted as aris- 
ing from the recoil the atoms undergo when they acquire 
mhk of momenta by scattering from the MS. 

The classical motion of the atoms in the y-direction 
associates a distance y with a time 

t = y/u, (25) 

where u = p y /M and M is the atomic mass. For a given 
it, the momentum space wave function evolves as 

^(p,t)=e- i ^ h ^(p) 

= -\/27rfi,y^ rj m exp[— im 2 u>kt]8{p — mhk), (26) 

m 

where e p — p 2 /2M is the kinetic energy of an atom hav- 
ing center of mass momentum p, and u>k = hk 2 /2M is a 
recoil frequency. In the coordinate representation 

#M)= r -4=e^/ h 4>(p,t) (27) 
■/-oo V 2nh 

one finds 

ip(x, t) = rj m exp[imkx — im 2 (j> t ], (28) 

m 

where the Talbot phase is defined by 

<k = Ukt. (29) 

Superposition of the different terms in Eq. ( p8| ) leads 
to a spatial modulation of the atomic density 

f{x,t) = \4,{x,t)\ 2 . (30) 

The interference terms in Eq. ( |30| ) are a direct mani- 
festation of matter-wave interference. One can see that, 
as a function of the time of flight t = y/u, the wave 
function (ESh undergoes oscillations on a time scale w^ 1 . 
As a consequence, the atomic spatial distribution (pG) 
contains quantum beats at frequencies (m 2 — n 2 )uik, for 
integral m, n. Such quantum beats have been predicted 
by Chebotayev et al. (1985) and observed by Chapman 
et al. (1995). 
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It follows from Eq.(p8|) that the atomic wave function 
coincides with the amplitude transmission function of the 
MS when 

t = t T = 2ir/u k . (31) 

At this time, atoms are found in the focal plane at 

y = L T = ut T = 2d 2 /\ dB , (32) 

and a self-image of the MS is produced. In general one 
finds that the atomic wave function is a periodic function 
of the Talbot phase <j) t having period 2tt, a periodic func- 
tion of the time t having period 27r/cjfc, and a periodic 
function of the distance y having period Lt- 

The self-imaging of a periodic structure is well known 
in classical optics as the Talbot effect. To describe this 
effect, one usually starts from the Fresnel-Kirchhoff equa- 
tion 

1 f°° 

ip(x) = / dx'rj(x') exp \ik dB {x - x') 2 /2y\ , 

VlAdBV J-oc 

(33) 

which is written here in the parabolic approximation. To 
establish an equivalence between Eqs.(|3^) and (p8|), one 
can substitute the Fourier expansion of the function ry(x') 
in Eq.(]3^), carry out the integration, and express y in 
terms of t. 

It is possible to derive a useful symmetry property for 
f(x,t) when the transmission function rj(x) is real, as it 
is for the MS. For real r](x), there is pure amplitude mod- 
ulation of the atomic wave function and r\ m — rf_ m . It 
then follows from Eqs. (]2^) and ( |30| ) that the atomic spa- 
tial distribution is invariant under inversion with respect 
to the plane y = Lt/2, i. e. 



(34) 



ip(x, <f> t ) = (1/d) / r)(x')Z(x - x', <f> t )dx', (35) 

J x — d 



where 



Z(x, 4>t) = exp [— im 2 (j)t + imkxj ■ 



(36) 



In the following discussion, we calculate ip{ x : 4>t 
2TTy/Lr) at fractions of the Talbot length, that is, for 



V = L T /n, 



or, equivalently, for 



<t> t = 27r/n, 



(37) 



(38) 



where n is a positive integer. 

When n = 2, one can show that the diffraction pattern 
is a self-image of the MS, shifted by half a period. For 
n = 2, 4>t = 7r, and 



Z{x,tt) = (l-e lkx )J2e 2tqkx - 

9 

Using the equality 

J2e tqa = 27r^(5(a-27rs), 



finds 



Z(x, it) — d 8 



odd 



X — S — 

2 



(39) 



(40) 



(41) 



For the integration range in Eq. (pq), only the s = 1 
term contributes when Eq. ( |4l| ) is substituted into Eq. 
(B5f), leading to 



ip(x,Tr) =T) ( x- - ) , 



(42) 



f{x,t)=f(x,t T -t) 



f (x,t)\= f (x,t)\ 



L T ~y ■ 



(34a) 



(34b) 



Thus, one need calculate / (x, t) in the range < y < 
Lt/2 to obtain the distribution for all y. 



B. Calculation of the atomic density profile 

We have seen that, for a sufficiently collimated atomic 
beam, self-imaging of a MS occurs at integral multiples 
of the Talbot length. To analyze the diffraction pattern 
for arbitrary y, it is convenient to use Eqs. (^4|) and ( psj ) 
to reexpress the atomic wave-function as the convolution 
(Winthrop and Worthington, 1965) 



i. e. at half-integral multiples of the Talbot length, there 
are self images of the MS shifted by half a period. 
For arbitrary n, it is convenient to write 



m = nq + r, 



(43) 



where < r < n — 1 and q and r are integers. It then 
follows that 



exp (— im 2 (j)t) — exp(— 2mm 2 /n) = exp (— i2i:r 2 /n) 



and 



where 



Z(x, 2-7T jn) — <f ^""^ a s S (x— s— j , 



(44) 



(45) 
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^ n— 1 

a s (ra) = — } exp [27rir(s 
n z — ' 

r=0 



-)/»] 



(45a) 



The atomic wave function ( |3£ 
n-i 



27r/n) = ^2 a s{n)rj 



is then given by 



s=0 



The meaning of this equation is clear. At distances 
y = Lt I n , the wave function consists of n self-images of 
the amplitude transmission function 77(2), having differ- 
ent amplitudes a s (n) (some of which might vanish) and 
spaced from one another by the distance d/n. For ap- 
propriately chosen r)(x) and n (see below), the atomic 
density (|30|) associated with the wave-function (^) is a 
periodic function of x having period 



d g = d/n, or d g — 2d/n. 



Thus, the Talbot-effect can be used to generate spatial 
modulation of an atomic beam having a period which is a 
fraction of the period of the MS. We refer to such profiles 
as "higher order atomic gratings." 



To simplify the expression for the coefficients a s (n), 
one can use an alternative approach for evaluating 
Z(x,cj)t) (Winthrop and Worthington, 1965). The sum 
in Eq. (pq) can be written in the form 



(46) 5; £/ 



1 

x exp 



dmdf 
if(q - m) 



4>t + imkx] . (48) 



Carrying out the integration over m, summation over q 
[using Eq. ([40|)1 and integration over /, one arrives at 



Z(x,(f> t ) 



■zr/~2 z r{x,<j)t), 

l(f>t ^ 



(47) where 



z r (x, 4>t) = exp i (kx + 2ttt) /4<p t 



At distances y — Lt / n one finds 



z r (x, 2ir/n) 



exp 



in (kx) /8ir 



exp \inkx (q + o) + in\ 1 
exp (inkxq) , 



for r = 2q+ 1, 
for r = 2q, 



(49) 



(49a) 



(50) 



where q is integer. Substituting this expression into 
Eq. ( p9| ) and summing over r, one arrives again at Eq. 
([45|), but with an alternative expression for a s (n): 



a s (n) 



1 



l + (-l) s e m7r/2 



2 /2r, 



(51) 



C. Higher-order gratings using the Talbot effect 



One can conclude from Eq. (|4^) that, owing to matter- 
wave interference, the transmission function 77(2;) im- 
printed on the atomic wave function by the MS can be 
copied n times in the plane y = Lt/ti, with each copy 
separated by d/n. This effect occurs for arbitrary trans- 
mission functions and can be used to generate higher or- 
der atomic gratings. 

Since the a s (n) appearing in Eq. ( f46| ) are not necessar- 
ily equal, the wave function ( |46"| ) is periodic with period 
d, but not necessarily with period d g < d. Moreover, it is 
possible for the different copies corresponding to differ- 
ent s to overlap. We refer to a pure, higher order atomic 
grating as one in which the different grating images do 
not overlap and for which d g < d. When the width fd of 
the slits in the MS is smaller than spacing d/n, different 
terms in the wave-function ( |46] ) do not overlap with one 
another and one finds for the atomic density (pC 



/(M)=£k HI 2 



77 1 x — s— 

n 



(52) 



Pure higher order atomic gratings are produced only 
if the nonvanishing |a s (n)| are equal. From Eq. (|5l| ) one 
sees that 



\a s (n)\ 



2 ( |cos(n7r/4)| , for even s 



n 1 I sin (n7r/4)| , for odd s 
Three different situations can be distinguished 

1. n = 2m+ 1 

2. n = 2(2m+ 1) 

3. n = 4m 



(53) 



(54) 



for integers m > 0. 

In the first case |a. s (n) | = -4=, independent of s. At 
distances y — Lt/3, y — Lt/5, . . ., pure, higher order 
atomic gratings having periods d g = d/3, d/5 ... are 
produced. 

In the second case |a s (?i)| = for s even and |» s (n)| = 

independent of s, for s odd. In the plane y — L^/2 

the atomic grating is shifted by a half-period d/2 from the 
initial grating (as found above); in the plane y — Xy/6 
only terms located at x = 4, 



^ in Eq. @ con- 



tribute to the sum (for < x < d), corresponding to an 
atomic grating having period d g = c?/3, that is shifted 
by a distance d/6 from the initial grating. In general in 
the focal plane y = Lt/[2 (2m + 1)], one finds an atomic 
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grating having period d g = dj (2m + 1) which is shifted 
by a distance d/[2 (2m + 1)] from the initial grating. 



IV. SHADOW EFFECT WITH 
MICRO FABRICATED STRUCTURES 
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FIG. 3. One period of the atomic density spatial modula- 
tion in the planes located at the distances y — Lt /m (m = 1, 
. . . 16). The case of the atomic beam modulation with micro- 
fabricated structure having transmission step-function ( p?j[ ) 
with relative width / = 0.16 is shown. While expected period 
d g of the beam self-image greater than slits' width a = df, 
one obtains higher-order spatial gratings. 



In the third case |a s (n)| = y ^, independent of s, 
for s even and |a s (n)| = for s odd. In the planes 
y = £t/4, y = Lt/8 . . ., one finds atomic gratings hav- 
ing periods d g = d/2, d/A . . .. 

The atomic density profile can no longer be written in 
the form (p2J) when 



In the previous section, it was assumed that the an- 
gular divergence 9b of the incident beam was less than 
d/LiT- If this inequality is not satisfied, the diffraction 
patterns associated with different velocity subgroups in 
the incident atomic beam result in a washing out of the 
overall diffraction pattern. For typical beam parameters, 
this condition restricts 9b to be less than 10~ 5 — 1CP 4 rad. 
The restriction on 9b is a limiting factor on the maxi- 
mum flux of the atomic beam. It is possible to avoid 
this restriction and increase the atomic flux if echo-like 
techniques are used. 

Using echo techniques that are analogous to those en- 
countered in the study of coherent transients, one can 
observe matter- wave interference in beams having a large 
angular divergence (Dubetsky et al., 1984). It turns out, 
however, that the dephasing and rephasing of the atomic 
gratings which occurs in such schemes does not depend in 
any critical manner on quantization of the atomic center- 
of-mass motion. In other words, the dephasing-rephasing 
mechanism is the same whether or not L ~ Lt (Talbot 
effect) or L <C Lt (classical limit). As such, it makes 
sense to consider the limit of classical scattering first, 
since the analysis is easier and a simple geometric in- 
terpretation can be given to the results (Dubetsky and 
Berman, 1994). Thus, we consider the limit L <C Lt in 
this section and defer a discussion of the case L ~ Lt 
until Sec. V. 

In this and the following section we consider the inter- 
action of an atomic beam with two MS, separated by a 
distance L. The angular divergence of the incident beam 
is sufficiently large to satisfy the inequality 



9 b > d/L. 



(57) 



fd > d g 



(55) 



In this limit, different components in the wave- function 
( |46| ) overlap and can interfere with one another in form- 
ing the atomic density. Even though the atomic distribu- 
tion function can still contain narrow peaks having a size 
of order of d g , the amplitudes of the peaks are not equal, 
and the period of the overall diffraction pattern reverts 
to the period d of the initial grating. 

These different regimes are illustrated in Fig. || plot- 
ted for an amplitude transmission function defined in the 
interval < x < d as 



r](x) 



1 . for < x < df 

0, for df < x < d 



(56) 



The Talbot effect enables one to create pure, higher-order 
atomic gratings having periods that are limited only by 
the slit widths in the MS. 



The first MS produces a sum of atomic gratings, one for 
each velocity subgroup in the initial atomic beam. Imme- 
diately following the MS, these gratings overlap and mir- 
ror the transmission function of the grating, but down- 
stream from the MS, they dephase relative to one an- 
other. As a result, the macroscopic atomic grating is 
washed out at a distance 



1 - d/9 b < L 



(58) 



from the MS. 

Although the macroscopic grating produced by the 
first MS washes out in a distance of order I L, it 
is possible for the second MS to lead to a restoration of 
the atomic gratings. For particles moving on classical 
trajectories, we refer to this process as a shadow effect 
since it can be interpreted completely by the "shadow" 
of the incident beam formed by the two MS (Chebotayev 
et al., 1985) (see Fig. |). 
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Microfabricated Atomic gratings 




™™ ' ■► y Gratings' periods 

FIG. 4. Shadow effect, geometric simulation. Atomic 
beam having large angular divergence (Ot = 120°) passes 
through two separated MSs having the same periods 
di — d,2 — d, duty cycles /i = 0.5 and fa = 0.05 and 
separated on the distance L. Grating carved by first MS is 
washed out on the distance £ ~ d. Other gratings between 
MSs (0 < y < L) are Active, they appeared because only 
trajectories for atoms passed through both MSs are shown. 
Gratings after the second MS (y > L) are real, they appear 
owing to echo effect. 

For a beam having large angular divergence (0b ~ 1), 
the initial grating produced by the first MS washes out 
in a distance comparable with the MS's period, in accor- 
dance with Eq. (|5^). After passing through the second 
MS, however, macroscopic gratings reappear in specific 
focal planes. A grating having the same period as the MS 
is focused in the focal plane y = 2L, while higher-order 
gratings having periods d/n (for integer n) are focused at 
other locations (to be determined below). The shadow 
effect can also be demonstrated easily using incoherent 
light (Chebotayev, 1986). 

Although the shadow effect occurs for classical parti- 
cles, it can be interpreted in terms of a dephasing and 
rephasing of atomic gratings. The relevant phases are 
the Doppler phases associated with various Fourier com- 
ponents of the atomic density, as discussed in the In- 
troduction. In such a picture, the final image on the 
screen depends on a cancellation of Doppler phases in 
the spatial regions — > L and L — > 2L, for example. In 
other words, the signal is sensitive to the relative Doppler 
phases in two spatial regions and is a measure of this 
relative phase. Insofar as interferometers are measures 
of relative phase, the echo-like rephasing of the atomic 
gratings can be viewed as a manifestation of atom in- 
terferometry. On the other hand, this rephasing is not 
related to the wave nature of matter. A shadow effect 



interferometer of this type was used by Batelaan et al. 
(1996) to measure the displacement of atomic gratings 
produced by rotation and by gravity. 

The same type of Doppler dephasing and rephasing 
that occurs using MS can also occur when atoms in- 
teract with two or more nearly resonant standing wave 
fields (Baklanov et al., 1976; Barger et. al., 1979; Dubet- 
sky, 1976; Chebotayev, 1978; Chebotayev et al., 1978a, 
LeGouet and Berman, 1979; Mossberg et al., 1979; Du- 
betsky and Semibalamut, 1982; Borde, 1989; Dubetsky 
and Berman, 1994). When standing wave optical fields 
are used for modulation of the atomic spatial distribu- 
tion, the atomic gratings often are monitored by applying 
a probe pulse in the focal planes that transfers the phase 
associated with an atomic state population to one asso- 
ciated with an atomic coherence. Atom interferometers 
of this type have been used for precision measurements 
of gravitational (Kasevich and Chu, 1991) and inertial 
(Riehle et al., 1991) phenomena [for a review, see (Mullcr 
et al., 1995)]. In these cases, external fields give rise to a 
displacement of the atomic gratings. 

A. Dephasing-rephasing processes using two 
spatially separated MS 

Before calculating the particles' distribution function, 
we derive some general properties of grating formation. 
In this subsection, it is convenient to make a Fourier 
decomposition of the atomic density profile in the x- 
direction. The propagation of each of the Fourier com- 
ponents is then treated separately. 

Consider the case when the two MS (1 and 2) have 
periods d\ and c?2 and are separated from one another 
by a distance L. A MS forms a periodic spatial distribu- 
tion (shadow) which is the same for all atomic velocity 
subgroups just after passing through the MS. The profile 
created by the first MS contains a sum of harmonics in 
the x— direction having spatial periods 

d mi = d/ \mx \ , (59) 

where mi is an integer. Immediately following the MS, 
the mith spatial harmonic varies as cos(m\kix), where 
k\ = 2n/di is the wave- number associated with the first 
MS. As the atoms move downstream from the first MS, 
the mith harmonic acquires a Doppler phase (|7a|) given 
by 

(pmi(t) = mihvt, (60) 

where v is the x-component of atomic velocity and t = 
y/u as before. For a time 

t d ~ l/{k x v) (61) 

the Doppler phases becomes large, 4> mi (t) > 1, and 
the macroscopic grating washes out on averaging over 
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v. Since v ~ 
distance £ (p£ 



u0b, the time td (|61|) corresponds to the 



The atoms pass through the second MS at time T = 
L/u(y = L). Downstream from the second MS, each spa- 
tial harmonic acquires an additional phase (i.e., a phase 
in addition to <fi mi (t), which, itself, continues to increase 
following the second MS) 



<t>m 3 (t) — m 2 k 2 v(t — T), 



(62) 



where 777,2 is another integer. Since the mask created by 
the second MS is superimposed on the shadow from the 
first MS, the resulting shadow consists of harmonics hav- 
ing wave numbers 



kh = Imifei + to 2 /c2| 



and velocity-dependent Doppler phases 



(63) 



<P (*) = 4> mx (t) + 0rn 2 (t) • 



(64) 



The two phases in the rhs of this equation can cancel one 
another at the so-called echo time t e defined by 



<M*e)=0, 



(65) 



corresponding to a focal plane y e = ut e . At such times, 
one produces an harmonic in the atomic density that is 
independent of v; as a consequence this grating survives 
any averaging over the velocity distribution in the inci- 
dent beam. From Eq. (65) one sees gratings are focused 
when 



t e /T = y e /L 



1 + (mi/m 2 )(fci/fc 2 ) 



(66) 



<D(t) 



T WT 



t 



FIG. 5. A dephasing-rephasing process responsible for the 
atomic grating focusing at the focal plane y — produced 
by two MSs having periods di and d,2 and associated wave 
numbers fei and fe, such as jki = £k2, where j and I are inte- 
grals. When atomic beam passes through first MS a number 
of harmonics having wave numbers m\k\ (harmonics m\k\) 
are induced. Mask, consisting of harmonics mikv, is superim- 
posed on the atomic spatial profile when beam passes through 
second MS. Dephasing of the relevant harmonic j(m — n)ki 
occurs between MSs. Second MS starts rephasing process for 
harmonic j(m — n)ki — £mk2 = —jnki which leads to the 
phase cancellation at the echo-point t = —T, where macro- 
scopic atomic grating is focused. Other harmonics contribut- 
ing to this grating must have wave numbers q-times larger, 
where q is integral. Their phase diagrams differ from plotted 
only by scaling the phase in q— times. 



[To satisfy Eq. (|65|), the integers mi and mi must have 
opposite signs (for the sake of illustration, we take ra\ as 
negative in the examples given below)]. 

We assume that the ratio of the MS's periods is ratio- 
nal, 



di/d 2 = k 2 /ki = j/£, 



(67) 



where j and £ are the smallest positive integers which can 
be used to satisfy this equation. In this case an infinite 
number of harmonics associated with pairs (7711,7712), 
having the same ratio contribute at a given focal 
plane. It follows from Eqs. (66, ^7|) that for 



mi 

7772 



(68) 



where m and 71 are positive integers having no common 
factors with m > n, gratings in the atomic spatial distri- 
bution appear at echo times 



tr 



n 



(69) 



The dephasing-rephasing process is illustrated in Fig. ||. or at focal planes located at 
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m 

Ve = —L. 

n 



From Eq. (|63|), one finds that harmonics having 

Imilfci 



m ^ 



(70) 



(71) 



are focused in this plane. 

For example, consider the limiting case in which d\ = 
d 2 = d (j — I — 1), analogous to the situation stud- 
ied by Dubetsky and Berman (1994). In the plane 
y = 2L (m/n = 2), all harmonics having m 2 /mi = —2 
(i.e. {mi, 1712} = { — 1, 2}; {— 2, 4}; {— 3, 6};etc) are fo- 
cused. As a result [see Eq. (Ba)], harmonics having 
kh = k, 2k, 3fc,etc, are focused in the plane y = 2L. 
The period of this atomic grating d g corresponds to 



the smallest value of k^, namely d„ 



27r/fc 



d. 

Similarly, in the plane y — 3L/2 (m/n = 3/2), all 
harmonics having 7712/7711 = —3 (i.e. {7774,7772} — 
{— 1, 3h {— 2, 6}; {—3, 9};etc) are focused. As a result [see 
Eq. (p3|)l , harmonics having kh = 2k, Ak, 6fc,etc, are fo- 
cused in the plane y — (2>/2)L. The period of this atomic 
grating is d g — 2ir/2k = d/2. For 777 = (77 + 1), one finds 
that atomic gratings having period d g — d/n are focused 
in the plane y = [(77 + l)/n]L. 

To treat the case of arbitrary, rational d\/d 2 = j/£, we 
set 

777,1 = — j{m — n)q, 7772 = imq, (72) 
where q is an integer, 

j = j/», e = t/n, (73) 

and [i is the largest common factor of j(m — n) and £m. 
Harmonics having wave numbers (|7l|) 



k h =njki\q\ 



(74) 



are focused in the plane y e = ^L. The minimum possible 
wave number 



k g = njkx 

determines the period of the focused grating 



2tt 



jn' 



(75) 



(76) 



One concludes that it is possible to create a higher 
order atomic grating, having a period that is jn— times 
smaller than that of the first MS, by passing an atomic 
beam having a large angular divergence through two 
MS. Although both the Talbot and shadow effects lead 
to higher order atomic gratings, there is a qualitative 
difference between the two cases. In the Talbot effect, 
the structure of the MS is copied 77 times in the image 
plane y — Lr/n, giving rise to a profile having period 



d g = di/n or 2d\/n, provided that d g > f\d\. The mini- 
mum period is determined by the slit width. In contrast, 
the period of the atomic grating produced by the shadow 
effect in the plane y e = is given by Eq. ([76]) and 
is not limited by the slit widths of the MS (although 
the contrast is determined by the slit widths). The pe- 
riod of the atomic grating is equal to d\/j in the focal 
plane y e — 2L and is compressed by a factor 77 in the 
plane y e = —L. This compression lies at the heart of the 
shadow effect's application to atomic lithography (Du- 
betsky and Berman, 1994). Atomic gratings having peri- 
ods smaller than those of both the MS and even smaller 
than the slit widths of the MS can be obtained. In this 
respect the shadow effect has yet an additional advan- 
tage over the Talbot-effect, where higher order grating 
production is not accompanied by compression. 

Before proceeding to calculate the atomic density dis- 
tribution, we should like to estimate the depth of focus of 
the various gratings. The distances between focal planes 
are comparable with the distance L between the MS. One 
can estimate the depth of focus l g from the requirement 
that the phase ( |64| ) be smaller than unity in the region 
of the focal plane. For mi given by (|72|) one finds 



4> (t) = qkgvSy/u, 



(77) 



where Sy = (y — —L) = (y — y e ) in the neighborhood of 
the focal plane at y = ~L. Setting q = 1, <j){t) ~ 1 and 
6y = £ g , and using Eq. ( |75| ) and the fact that 8 b ~ v/u, 
one obtains 



— <i, 
nj 



(78) 



where t is given by ( p8[ ) with d = d\. Since L I has 
been assumed, it is possible to separate the various grat- 
ings. The sharpening of depth of focus of the higher order 
gratings predicted by Eq. ( |78|) is in qualitative agreement 
with the results shown in Fig. |]. 



B. Particles' distribution profile 

We now turn our attention to a calculation of the 
atomic density profile. In contrast to the Talbot effect, 
it is not possible to find self-imaging using the shadow 
effect, since the scattering coefficients for the different 
spatial harmonics are not the same. Let the transmission 
functions for the two MS be denoted by Xs{x) (s = 1 or 
2) [xs{x) is a transmission function for atomic density, 
while n s (x) is a transmission function for atomic state 
amplitudes - for MS having transmission of either 1 or 
0, these functions are identical]. Atoms are scattered by 
the MS in the planes y = and y = L, or, equivalently, 
at times T\ = and T2 = T = L/u. Calculations are 
carried out using t = as a variable. In some sense, 
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this corresponds to working in the atomic rest frame. As 
a result of scattering, the atomic density is modified as 



f(x,v,T+)=Xs{x)f(x,v,T s ), 



(79) 



where Tj^ are times just after or before a scattering event. 
Following the scattering event, the distribution evolves as 



f(x, v, t) = f[x - v(t - T s ), v, T s 4 



(80) 



We assume that, for t < 0, the atoms are distributed 
homogeneously in the transverse direction, i.e. 



where 4>(t), as defined by Eq. ( |64|), is also a function 
of m and n. Owing to condition (|57|), for t ~ T ~ L/u 
and v ~ u9 0l the phase factor in the brackets of Eq. (p3) 
is large, of order L9 /d s 3> 1. On averaging over veloc- 
ities, one finds a nonvanishing contribution only at the 
particular focal planes or echo-times given by Eq. (65). 

Retaining contributions from only those rrii given by 
([72|) corresponding to the various focal planes, one finds 
from Eq. (H ) that the atomic density in the focal planes 
is given by 



f(x,v,t)\ t 



<o 



1. 



(81) 



f(x,t e 



q 



(2) 



(85) 



The assumption of a homogeneous velocity distribution 
is consistent with a beam having angular divergence (57), 
since in this limit, the transverse velocity distribution is 
approximately constant over the range d s /L. The spatial 
distribution of the atomic density for t > T is given by 



where k g is the wave number of the focused atomic grat- 
ing given in Eq. (f75|) . Thus, owing to the shadow effect, 
at the echo-time (69) an atomic grating is focused having 
period 



f(x,t) = (xi(x - vt)x2[x - v(t - T)]) 



(82) 



d g = 2n/k g = di/(jn). 



(86) 



where (. . .) represents an average over velocities. Ex- 
panding Xs (x) in a Fourier series 



(s) imk s x 



(83) 



where k s — 2n/d s is the "wave-number" of structure s, 
one finds 

fi x > l ) = X! Xm\Xm 2 exp[i(fcimi + k 2 m 2 )x\ (exp [-#(*)]) 

mi, 7712 

(84) 

fdi f d 2 d ^d>r.„ „ . . 

f(x,t e ) -- 



Note that the period of this grating is nj times smaller 
than the period of the first microfabricated structure. 

To reexpress the density function at the focal planes in 
terms of the transmission functions, we write the Fourier 
harmonic amplitude Xs as 



X 



U) 



^rXj(x)exp(-ikjSx), 



(87) 



from which one finds the atomic density at the focal 
planes given by 



' ' 2 ^ exp {iq [k g x + j (m - n) k x x x - m£k 2 x 2 ] } Xi ( x i) X2 (x 2 ) 
9 



d\d 2 



The sum over q leads to an infinite set of 5— functions 



9 



a iqa 



27r^(5(a-27rs') 



which allows one to carry out the integration over xi, for example. Only the values 



x\ 



di 



j(m — n) 

in the range [0, di) contribute, resulting in the inequality 



-x 2 x 
di dn 



X « 2 ^ I X 7>X2 , — / 

mi — < s < mi h Jym - 

d g d 2 d g d 2 

for integer s' . One finds that a finite number of terms contribute, having 

m£-r- — 



(89) 



(90) 



(91) 



(92) 



where s = 0, 1, ... j(m — n) — 1. In Eq. (92) and for the remainder of this chapter, a notation is adopted, in which 
we set A = [A]j + {A} F , where [A]j and {A} F are the integral and fractional parts of A, respectively. Using Eqs. 
@ and @, we obtain 
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^ j(m-n)-l d2 
f(x,t e ) = T r V / ~TX2 (x 2 ) XI 

j(m-n) J d 2 



di 



j(m — n) 



ml— - — 

d 2 d 



a J f 



(93) 



Consider in detail the case when the micro-fabricated structures have duty cycles (ratio of slit openings to periods) 



fj and 

for {#1 

Introducing dimensionless variables 

w = x/dg, z = x 2 /{d 2 f 2 ) 
and taking into account that the argument of function xi m Eq. (B3h is positive, one obtains 



(94) 



(95) 



/(Me) 



7 (to — n) / — 

J y ' s=0 



where 



h s (w) — / dzO [(3 — (s + {az — w} F )] , 
Jo 



(96) 



(96a) 



a = m£f 2 , P = j(m-n) fi, 

1 / ov x ^ 

and 9(x) = i ' . _ is the Heaviside step function. It is sufficient to consider only the range 

1 0, for x < y 

< w < 1. 



For < s < \J3]j — 1, the integrand in Eq. (96a) is equal to unity. Therefore, 

h 



f(x,t e ) = 



j (to — n) 



and one needs to evaluate the expression (96a) only for s 



VhW - / " 1 i 'i r 

10 



dz0({(3} F - {az - w} F ) 



(96b) 



(97) 



(98) 



(99) 



The first term in the Eq. (|98|) brackets is independent of w = consequently, it is only the second term which 
corresponds to the atomic gratings. 

A method for evaluating the integral (K3) is given in the Appendix. Using this method one finds 



h mi (w) = [{P} F [a] z + S{w)]/c 



where the function S (w) is given by 

S(W): 



{a} F -w, for < w < 1 - {f3} F 

{(3} F + {a} F - 1, for 1 - {{3} F <w< {a} F 

{/3} F + w-l, for {a} F < w < 1 + {a} F - {13} 1 

{a} F , for 1 + {a} F - {13} F < w < 1 



(100) 



if {[3} F >max({a} F ,l-{a} F ), (101) 



S{w) = < 



{a} F — w, for < w < {a} F 

0, for {a} F <w<l- {[3} F 

{(3} F +w-l, for 1 - {(3} F <w<l + {a} F - {(3} 1 

, {a} F , for 1 + {a} F - {(3} F <w<l 



> if {a} F < {P} F <l-{a} 



F ■■ 



(101a) 
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S(w) 



{a} F - w, 
{13} F + {a} F 
{p} F + w-l, 



- 1. 



for < w < {a} F - {P} F 
for {a} F - {[3} F <w<l 
for 1 - {/3} F < w < {a} F 
for {a} F < w < 1 



if 1 - {a} F < {/5} F < {a} 



F ' 



(101b) 



{(3} F , for < w < {a} F - {/3} F ( 

{/?} F +«;-!, for 1 - {/3} F <w<l ' 



Substituting Eq. (100) in the rhs of Eq. ( pq ) one finds 
that the atomic beam profile at a given focal plane is 
equal to 



f(x 1 t e ) = [a[p} I + {p} F [a]! 

+ S (w — x j d g )] I [m(m — n)£j] 



(102) 



C. Main features 



All dependencies in ( |l0l| ) coincide if 
when 

S{w)-- 

In this case the grating amplitude 

A f\pO? ^e)max f{p*i 



(103) 



(104) 



(105) 



for given m,£,j,n, achieves a maximum value A = 
(2m£j(m — n)) -1 . To maximize this quantity for a given 
grating period, one has to choose I = j = 1 and m = n+1 
(£ = j = 1), which corresponds to the focal planes 
y = 2L (n = l),y= 3 / 2 L (n = 2),y= 4 / 3 L (n = 3) . . ., 
where gratings having periods d g = d\/n are focused. To 
satisfy condition (|103[), one can choose 



for which 



A = 1/2, / a = l/[2(n + l)], 



A = l/[2(n + l)]. 



(106) 



(107) 



The constant background term of f(x, t e ) [first two terms 
in the numerator of Eq. (102)] vanishes since 



3 



X 

a 

o 

3 

b 



Efl 
C 

o 

T3 



O 
< 



0.25 



| x=0 



y e =2L 



x=dj 



0.25^ 



i x=0 



y e = 3 / 2 L 



x=d! 



0.25 



i x=0 



y e =/ 3 L 



x=d[ 



0.25. 



| x=0 



y e =/ 4 L 



x=d[ 



0.25^ 



| x=0 



y e =/ 5 L 



x=d[ 



X 



FIG. 6. Gratings' profiles produced by two MSs having 
equal periods (di = e^) and duty cycles /i = J, = tk 
at the different focal planes. Chosen case is optimum for 
5— order grating focused at the plane y e = |L. This grat- 
ing has amplitude A = ^ and twice large background term. 
The amplitude of this grating is not less than amplitudes of 
the lower-order gratings focused at the other planes for given 
duty cycles. 

To achieve this maximum signal, one must use slits 
in the second microfabricated structure whose width is 
smaller than the atomic grating period {fidi = dx/2(n + 
1) < d g = dx/n). Using the shadow-effect technique, one 
can al so o bserve atomic gratings having the same ampli- 
tude (107) whose period is smaller than the slit width, 
but some background term appears for these gratings. 
Indeed, if = 2 (n+i) ' ^ — j — l,m = n + I (leading to 
{a} F = 1/2, (3 = 1/2), for positive integers q < n, then 



Eq. (103) still holds and provides the grating amplitude 
(107), but the background term is q times larger than the 
grating amplitude. 

For illustration we plot in Fig. ^| the grating pro- 
files at different focal planes y e = ^-j^-L (n — 1, . . .5) 
for MS having the same periods d\ — d 2 , and duty cy- 



(108) cles fx 



_ 5 



such that the parameters a and 
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P are given by a = y|(n + 1), (3 — |. This case corre- 
sponds to the 5th order grating at the plane y e = |L 
having maximum amplitude. When the amplitude of the 
5th order grating is optimized and / 2 ~ 0.5, the am- 
plitudes of the gratings that are focused in the planes 
y e = ^-^-L (n = 1, . . .4) are less than or equal to the 
amplitude of the 5th order grating. This feature is seen 
in Fig. | 

The geometric simulation introduced above allows one 
to obtain the positions of the atomic gratings. It can 
also be used to provide some quantita tive results. For 
example, one finds from Eqs. (101, 102 ) that the shadow 
effect disappears at the focal plane y = 2L if both MS 
have the same period (d\ = d,2 = d) and duty cycles, 
fi = $2 = 0.5. In th is case, m — 2, n = j = £ — 1, 
and, from Eq. fl96b| ), one finds that a = 1, /3 = 0.5. 
From Eq. ( 101a ) one finds that S (w) = 0; there is no 
atomic grating. The reason for the absence of the grating 
under these conditions is evident from Fig. 0. 

Microfabricated 
structures 



Atomic 
beam 





B 




: — — - 


0.5(1 


0.5d 




\ 



-0(x,2L) 



y=0 y=L 
FIG. 7. To explain an absence of shadow-effect at the focal 
plane y — 2L for half-open MSs having the same periods. An 
illumination of an arbitrary point O at this focal plane from 
given slit AB of the second MS is determined by the number 
of particles moving into the point O inside angle AOB, i. 
e. it is proportional to the length of the bold part of A'B' , 
given by s = \A'B'\ - f . Since \A'B'\ = 2 \AB\ = d, s is al- 
ways equal to 4 independently on the point O x— coordinate. 
Consequently, any variation of the particles' distribution at 
the focal plane y — 2L disappears. 



The geometric picture can also be used to explain 
the absence of background terms at the focal planes 
y = 2ii L produced by MS having duty cycles (106) 
and equal periods {d\ — g? 2 = d). One can see from Eq. 
(102, 104) that the background disappears because there 
are no particles at the points x q — ^ (q + |) , where q 
is integral. A geometric interpretation of this result is 
presented in Fig. ||. 
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FIG. 8. To prove that two MS having equal periods 
(da = d'2 = d) and duty cycles /i = | and /2 = 2(n+i) ( case 
n = 3 is shown) produce a background-free atomic grating 
at the focal plane y = ^=-L, one can notice that points 
x q = d g [q + |) are not achievable for particles (d g = ^ is a 
grating period). It follows from the fact that an angle, built 
from any point x q and arbitrary slit of the second MS, meets 
the closed part of the first MS. 



V. TALBOT-LAU EFFECT 

When the spatial separation of the MS is increased to 
the point where 



L - Ln 



2d 2 /X 



dB, 



(109) 



it is no longer possible to neglect quantization of the 
atomic center-of-mass in calculating the transverse mo- 
tion of the atoms. Just as in the Talbot effect, the recoil 
an atom undergoes on scattering from a grating must 
be taken into account. It turns out, however, that the 
Doppler dephasing and rephasing encountered in analyz- 
ing the problem of classical scattering by MS still can 
be given a classical interpretation when L ~ Lt, pro- 
vided that the angular divergence of the incident beam 
is sufficiently large, 0b 3> d/Lx = ^ds/d = 6d- In other 
words, even though we must account for quantization 
of the atoms' center-of-mass motion, effects related to 
Doppler dephasing (which are automatically included in 
a quantized motion approach) are unchanged from the 
classical case. 

We have already alluded to this result in the Introduc- 
tion. Recall that matter-wave interference results from 
the overlap on the screen of atomic wave functions asso- 
ciated with states having center-of-mass momenta p and 
p + nhk. The relative dephasing between these states 
M) contains a Doppler part (J7a) and a quantum part 



1G 



( ffh| ) . Since the washing out and restoration of the macro- 
scopic atomic gratings is connected with an averaging 
over atomic velocities v, one expects that the Doppler 
part, proportional to v, is responsible for the dephasing- 
rephasing effect. This contribution is actually classical 
in nature (i. e. it does not vanish in the limit Ti — > 0) 
and enters the calculations whether or not the quantum 
contribution to the phase has to be considered. As a con- 
sequence, the dephasing-rephasing process is the same for 
the classical shadow effect and the quantum Talbot-Lau 
effect. It turns out, however, that for separations of the 
MS equal to a rational multiple of the Talbot length, 
the Talbot effect can actually result in a decrease of the 
period of the atomic gratings from those periods which 
would result from the classical shadow effect result. The 
decrease in period occurs for MS consisting of open slits 
and opaque strips; it would not occur, for example, in 
resonant standing wave fields. 

Since the Doppler dephasing determines the position of 
the focal planes and the period of the atomic grating, one 
can carry over the results Eqs. (|9[ |7(], f76|) obtained in 
Sec. IV for the shadow effect. In this section, we are in- 
terested in the variation of the atomic gratings in a given 
focal plane as a function of the separation of the MS. In 
other words, we look for those separations L for which 
the Talbot effect significantly modifies the gratings that 
would have been produced by the shadow effect alone. 
This is analogous to photon echo studies of atomic re- 
laxation in which the echo amplitude is monitored as a 
function of the separation between the excitation pulses. 

It should be noted that the Talbot-Lau effect has been 
studied using light by Clauser and Reinsch (1992) for the 
parameters 



di 



= 3; y e = 3L, 



(110) 



corresponding to {m, n, j, £,j,£}={3, 1, 3, 1,1,1/3} 
[recall that j = j/fJ,;£ — £/p., where /i is the largest 
common factor of j(m — n) and ml] and a grating pe- 
riod dg — d\. The atomic Talbot-Lau effect was demon- 
strated by Clauser and Li (1994) using K atoms for the 
parameters 



di 9 

-T- = 2, y e 

d 2 



2L, 



(111) 



{to, n, j, £,],£} ={2, 1, 2, 1,1,1/2},^ = ^. A the- 
oretical study of the atomic Tal bot-L au effect was also 
carried out for the parameters ( 111 ) by Carnal et al. 
(1995). The conclusions as to the period and location 
of the atomic gratings follows from purely classical con- 
siderations in this case; there is no need to invoke argu- 
ments related to the wave nature of matter (Clauser and 
Reinsch, 1992; Clauser and Li, 1994; Carnal et al., 1995). 



A. Grating formation 

The geometry is the same as that considered for the 
classical shadow effect, except that L is no longer re- 
stricted to be less than Lt- Again, it is convenient to 
work in the atomic rest frame defined by t — y/u. As 
discussed above, it is necessary to quantize the atomic 
motion only in the x-direction. The atoms undergo scat- 
tering at the MS at times T s = L s /u (T x =0, T 2 = T = 
L/u). For thin gratings, the atomic wave function i/j(x, t) 
undergoes jumps at the MS given by 



ip(x,T+) = Va (x)tp(x,T a ), 



(112) 



where i/)(x,Tg) is the wave function on either side of 
grating s, 7] s (x) is the amplitude transmission function 
of grating s, and 



Xs{x) = \r)s{x)\* 



(113) 



is the transmission function of grating s [for MS consist- 
ing of a series of slits, r) s (x) = Xs(x)]- 

To characterize the atomic beam using a quantized 
center-of-mass description, one can use the Wigner dis- 
tribution function defined by 



f(x,p,t) 



d.f 
2nh 



e^p(—ipx/h)ip(x - 



-,t)ij)*{x- -,t). 

(114) 



For scattering at a MS, one finds 



f(x,P,T+) 



dxdp' 
2nh 



exp [-i (p -p')x/h]r) s ( x + - ] //.; 



^)f{x,p',T; 



(115) 



When T] s (x) is a periodic function of x, one can write Eq. (115) as 



f{x,P,T+)=^2cxp(im s k s x)rj < n s ] rff f 



fik i 
x,p - —(n s +n s ),T s 



(116) 



where 



(116a) 



and 
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(116b) 



is a Fourier component of rj s (x), having period d s and wave number k s — 2it /d s . For times other than T a , the Wigner 
distribution function evolves freely as 



f(x,p,t) = f(x - v(t - T s ),p,T+), 

where v = p/M . 

Applying Eqs. (11161, |ll7j) one obtains the atomic distribution function for times t > T (y > L) to be 



(117) 



f(x, P ,t)= ^ 



„(2) 
ln 2 



(2) 



exp ^ im\k\ 
hk 2 



> - r[t-T)- [v-^(n 2 + n' 2 ))T 



+ im 2 k 2 [x — v{t — T)]} f |x 
V - t| («i + ™i) + fc 2 («2 + n' 2 )] 



v (n 2 + n' ) 



T-v(t-T), 



(118) 



where f(x,p) is the Wigner distribution function of the incoming atomic beam. The atomic spatial distribution is 
given by 



f{x,t) = I dpf {x,p,t) 



(119) 



which can be obtained from Eq. (118) as 



f( X ,t)= J2 ^ 



„(2) 



(2) 



dpf < x 



v + — (fci {m + rii) + k 2 (n 2 + n 2 )) 



(t-T) 



hki , 

V+ 2M (m+ni) 

Mi , ' 

v H rii + rii 



T,pj exp |z (mi/ci + m 2 k 2 ) 
[mikit + m 2 k 2 (t - T)\ 



hk 2 

X ~2M ( n2 + fl2 > 



(t-T) 



(120) 



In this expression terms having (mifci + m 2 k 2 ) ^ contribute to the atomic gratings. Owing to the assumption of 
an incident beam having large angular divergence Ob ~ v/u 3> d/L, the Doppler phases associated with these terms 
oscillate rapidly as a function of p, except in the echo focal planes. As a consequence, the positions and periods of 
the atomic gratings are the same as those in the classical shadow effect [see Eqs. (^, |7(], Uw]. In the remainder of 
this section, we calculate the atomic density in the focal planes y e or, equivalently, at times t e = y e /u given in Eqs. 

It is possible to simplify Eq. (120) if we assume that the angular divergence Ob of the incident beam is less than 
Oq = D/L, such that a freely propagating beam would undergo negligible diffraction over a distance of order L. For 
p ~ MuOb -C MD/T, one can neglect the dependence on m and n i of the distribution function appearing in Eq. 
(120). Then the sum over rii can be carried out using the formula 



E 

n 

-ivex/2 



-inot (s) 
° 'In 



'In — v 



a 

2k s 



Va[ x + 



2k, 



where 



r ds dx 

[F(x)] v = I -e-^F(x) 



(121) 



(122) 



is a Fourier component of the function F (x) . As a result one finds that the atomic density in the echo focal planes 
is given by 



f(x,t e ) = f(x)J2e iqkgX X ( ll 



j(m-n)q 



rfclx- qd 2 



4> T (m,n) 
2tt 



r] 2 [x + qd 2 



4> T {m,n) 
2ir 



(123) 



m£q 
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where Xs is a Fourier component of the transmission 



function xi {%) , 



(j> T {m,n) = 



j 2 (m — n) 
I 



(124) 



is a Talbot phase associated with a specific focal plane, 
and 



f(x)= / dpf (x,p) 



(125) 



is the initial spatial distribution in the atomic beam. 
Since the beam diameter is much larger than the period 
of the gratings, 



£> > d a 



(126) 



that, as defined by Eq.(127), there is a different Tal- 
bot length associated with the signal for different focal 
planes, y e = (m/n)L. We wish to examine the signal in 
a given focal plane as a function of the separation L of 
the MS or, equivalently, as a function of 4>t (tii, n) . When 
L = LT{m,n) \<j>r(m,n) = 2tt] the atomic density Jl23| ) 
is the same as that of the shadow effect (|S5|). 

Using arguments similar to those leading to Eqs. (j3j|), 
one can prove that, for pure amplitude modulation of 
the wave functions, i. e. for real amplitude transmission 
functions rjj(x) = r]*(x), the dependence of the parti- 
cles' distribution on the Talbot phase is symmetric with 
respect to the point </>y(m, n) = tt [L = Lxim, n)/2], 



(129) 



one can neglect the variation of / (x) and set / (x) = 1 
in Eq. (Of). 



The distribution function (123) is identical with the 
shadow effect result (§3), except for the presence of the 
Talbot phases. The main features of the dependence of 
the atomic density on the Talbot phase in the Talbot- 
Lau effect are the same as those for t he Talbot effect 
considered in Sec. III. The density (|l23| ) is an oscillating 
function of <^r(m, n) having period 2ir. If <j>T(m, n) is in- 
creased by 27r, or, equivalently, if the separation between 
the MS is increased by L<r(m, n), the density distribution 
in the corresponding echo plane is unchanged. The Tal- 
bot distance associated with a given focal plane is defined 
here as 



Lt(jti, n 



_ 2d 2 



XdB j 2 (m - n) 



(127) 



In terms of Lt(w, n), the Talbot phase (124) is equal to 



(^x(m, n) = 2 r K[L/Lx(rn, n)]. 



(128) 



In our notation, the Talbot phase is a function of L, 
while the Talbot distance is independent of L. Note 



f(x,t e )\ L = f(x,t e )\ LT _ 1 



(130) 



The question arises as to what values of <^r(m, n) lead 
to especially interesting results, i. e. atomic gratings 
that differ significantly from the gratings that would be 
produced by the shadow effect. We have found that the 
atomic gratings are significantly modified by the Talbot 
effect when the Talbot phase is a rational multiple of 27r, 



4>t(iti, n) = 2ir 



ttit 
tit 



(131) 



where and are positive integers having no com- 
mon factors. We proceed to analyze the atomic density 
function in the focal planes for separations of the MS cor- 



responding to Eq. (131), that is, for L = LT^^n)" 1 ^. 

For Talbot phases given by Eq. (131), the sum in Eq. 
(123) can be divided into ut — I independent sums having 



riTq + T 



(132) 



where < r < tit- For Talbot phases given by Eq. ( 131 ) 
any d ependence on q' disappears in the last factor of Eq 
(123), allowing one to rewrite Eq. (122) as 



TtT-l 



f(x,t e )= J2 E 



r— q' 



J 



dx\dx2 
d\d2 



exp {i (q'riT + r) [k g x + j (m — n) k\X\ — mifea^] } 



x Xi (zi) V2IX2- r d 2 7? 2 [x 2 + r d 2 , 

n T J \ n T 



(133) 



to sum over n' using Eq. (|89|), and to integrate over x%. The calculations are similar to those used to obtain Eq. (j9£ 
from Eq. (p8|), and one can obtain 



f(x,t e ) 



1 



jurim — n) 



n T -l j(m—n)n T -l rd 

E E 



dx 2 f 2nir 

— exp <^ 

«2 L n T 



s + 



tit ( m£ 



~x 2 x 



xr; 2 ( x 2 - r— d 2 ) tj 2 ( x 2 + r— d 2 ) xi , . 



di 



I I X2 x 
n T mi- — 

d 2 d„ 



(134) 
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B. Higher-order gratings using the Talbot-Lau 
technique 



Equation (134) is the basic result of this section. It 
gives the atomic density function in the focal plane for 
separation of the MS that corresponds to Talbot phases 
which are rational multiples of 2-7T. For specified trans- 
mission functions, it can be evaluated numerically in focal 
planes defined by t e = y e /u — (m/n)T = (m/n)L/u for 
arbitrary (j/£) = dx/d-2, (recall that j = — l/n, 

where /i is the largest common factor of j(m — n) and 
Im), rriT and riy [L = -Lt(w, n)^ 21 ]. In this subsec- 
tion, we are interested primarily in showing that, owing 
to the Talbot effect, periodic atomic density gratings can 
be produced whose periods d T are smaller than the cor- 
responding periods d g which would have been produced 
by the shadow effect. 

Th e first thing to note is that the function xi m Eq. 
( |l34| ), considered as a function of x, is periodic with pe- 
riod 



dg/riT, 



(135) 



nj— times smaller than the period d g of the shadow- 
effect grating. Unfortunately, this does not guarantee 
that f (x,t e ) is periodic with period d^, owing to the 
exponential term in Eq. (134). Under the transforma- 
tion x — > x + dT, the exponential term is multiplied by 
the phase-factor 



exp (2irir / tit) 



(136) 



which is a function of the summation index r. If the sum- 



mation over r in Eq. ( 134 ) somehow was restricted to 
r = 0, the atomic grating would have period dr- Re- 
stricting the summation to r = can be accomplished by 
choosing the amplitude transmission function such that 

the product 772 {%2 — ^f^^j V2 ( x 2 + 777^2) is nonvan- 
ishing only for r = 0. To simplify the discussion, we have 
taken m-r = 1 [L = LT{m,n)/riT]- 

For MS consisting of slits and opaque strips, both the 
amplitude transmission functions r\j (x) and transmission 
functions Xj( x ) are equal to the Heaviside step-function 



%O0 = Xj(x) 



a 3 ) F 



(137) 



The functions 



where fj is the duty cycle of MS j. 
1 23- ± ^ \ shown in Fig. [| represent the profile of 

do 



the second MS displaced by ± — «2 



{x 2 /d 2 +a} F 




Xo/cU 



f 2 +l-{a) F 



{x 2 /d 2 - a} F 



l-{a} 




w; ^h r^ x 2 /d 2 

FIG. 9. Plots of the functions ± 0} . In the shaded ar- 
eas only transmissions 4*2(3:2 ± acfo) 7^ 0. When slits relative 
width f2 < min ({a} , 1 — {a}) there are no shaded areas in 
the vicinity of the point X2 = 0. If also ft < 1 — 2{a}| , 
then shaded areas have no common points and product 
02 (£2 + ad2)4>2(x2 — ad2) = for any X2 



In the 

112 \ X2 - - 



range 
-d 2 ) rfc (X2 



< 



-d 2 



< 



the 



(/ 2 -{f 



product 

^) ) 

™ T J fJ 

8 ^/2 — I ^ — j J , which represents a product of 

profiles of the second MS displaced by ±^-d2, vanishes 
for r ^ provided that f 2 is sufficiently small and pro- 
vided that ^- 7^ i (if — I , the gratings are displaced 
by 4f and overlap for any ^2)- 



If 



f 2 < min ( — , 1 

tit tit 



(138) 



the only regions where 772(2:2 ± ^7(^2) does not vanish are 



02 



1 , /2+1 

TlT TIT 



and — £ 

d 2 



— J2 + — 

TlT TIT 



(139) 



respectively. These two intervals have no common re- 
gions if h + 1 - — < — or f 2 + — < 1 - — , i. e. 

o tit — tit J tit 

if 



1 



TlT 



(140) 



Inequality (14C) must hold for all r 7^ to guarantee that 
the atomic grating has period = d g /nr- Clearly, in- 
equality (14C) does not hold for r — ^ when tit is even. 
While this does not preclude the possibility of higher or- 
der gratings for tit even, it does suggest that we consider 
the cases of even and odd tit separately. 
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1. tit odd 

In this case, we write 

n T = 2n' + 1, 



(141) 



where n' is a positive integer or zero. For the summation 
range < r < — 1 in Eq. (134), the minimum value 
of the rhs of both inequalities (138) and (140) is l/rvr, 
which occurs for r = 1 or r = 2n' in ( |138| ) and r = n' or 
r = n' + 1 in (Pol). Thus, provided that 



1 



2ri + 1 



(142) 



one can produce atomic gratings having period dx = 
d g /(2n' + 1) in the focal plane y = for separations 
between the MS equal to L = ^n'+i) or > equivalently, for 
a Talbot phase (124) equal to t"ttt- Under these con- 



ditions one omits terms having r ^ in Eq. (134) and 
finds 



f(x,t e ) 



h 



j(m — n)riT 



[0% 



(143) 



h s {w)= dz9(/3' - (s + {a'z-w} F )), (143a) 
Jo 



a' = n T a = ml(2n' + 1) f 2 , 
0' = n T 0=J(m-n)(2n' + 



(143b) 



where a and (3 are given by (96b), and dimensionless 
variables 



x x 2 
d T d 2 j2 



have been introduced. Note that the ratio 

m f 2 i m f 2 d 2 



a/ (3 = a 1 10 



m-nfxj m-n f x d x 



(144) 



(145) 



depends on the focal plane and ratio of slit widths. In a 
manner similar to arriving at Eq. (102), one can obtain 



f(x,t) = [a , {p'} I + {f3 , } F [a'} I + S(w)} 



m(m — n)£j (2n' + 1)' 



(146) 



where S (w) is given by Eqs. (101 -101c) with the replace- 
ments a — ► a', [3 —> (3' . 




L=L T /5 






L=L T /4 



L=L T /7 




L=L T /8 
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L=L T /10 * L =L T /1 1 I L=L7 n 




rj o 



i I o 



x/d„ 



FIG. 10. Atomic spatial distribution created by two mi- 
crofabricated structures, having the same spacings di = d%, 
relative widths /i = |> /2 = \ and separated on the 
distance L. At the focal plane y = |L (m = 3, n = 2) , 
where owing to the shadow effect one expects grating with 
period d g = 4r, the wave matter interference leads to 
the higher order gratings if one chooses distance between 
fields coinciding with the integer fraction of the Talbot 



distance 



( L = Lt 



-, n T 



16 ] , which corre- 



sponds to the Talbot phase 4>t = For sufficiently 

small relative width of the second microfabricated structure, 
f'2 < -^p (/2 < ^7 J , one observes gratings having spacing 

dr = — = t^- I dr = — - = — 1 for distance between 

TIT Ztit \ tit ffp I 

fields equal to the odd (even) fraction part of the Talbot 
distance, i. e. for odd (even) tit- Exceptions here are the 
cases when tit is multiple 4 [L = L = L = -j^j , 
where for f2 < — wave-matter interference leads to the en- 
tire washing out of the gratings. 



The amplitude of the grating (146) is maximum when 
to = n + 1 , j = 1 = 1, and 



W} F = {13'} F = -, 



for which 



_ 2 gl + l 



2q 2 + l 



2(2n' + l)' J 2(n + l)(2n' + l) : 



(147) 



(148) 



where < q\ < 2n'. < q 2 < n are integers. Under 
these conditions one finds 



f(x,t e 



qiQ2 + \ (qi + 92) + 


f xn(2n'+l)\ x 







(n + 1) (2n' + 1)' 
This grating has amplitude 



(149) 
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A = \/ 2 (n + 1) (2n' + 1)' 



and a background term whose amplitude is (2q\q 2 + qi + 
52) times larger than A. Talbot-Lau gratings for different 
values of the Talbot phase are shown in Fig. HQ. 



2. rir even 

The atomic density patterns in Fig. |l^ have been 
drawn for both even and odd values of nr. From this 
figure one sees that qualitatively new features appear for 
even values of nr- When the Talbot phase (j>r(rn,n) — 
7r, 7r/3, 7r/5 (tit — 2, 6, 10) gratings having period 
d g , dg/3, dg/5, or 



(150) are focused. When the Talbot phase <pr(m,n) = 
7r/2, 7t/4, 7r/6 (nx = 4, 8, 12), the gratings are washed 
out entirely. To explain these results one needs to return 
to the general expression (134). If 

n T = 2ri, (152) 

where n' is a positive integer, one divides the sum over r 
in Eq. (134) into two parts having 

r = n'q + r', (153) 



d' T = 2d g /riT = dg/{nr/2) 
1 



(151) 



2jn l (m — n) 
■' v ' r'=0 s=0 



1' — 1 2j{m—n)n' — l 



with q = or 1 with r' restricted to the range < r' < 
n' — 1. In the second part (q = 1) one shifts the inte- 
gration variable from x 2 to (X2 + 5^2) j which leads to 
the same factors 772 , , Xi in both the q = and q = 1 
terms. In this manner, one arrives at the expression 



dx 2 
~d2~ 



x exp 



■Kir 



2n! ml 



x % [X2+r — d 2 j X i 



X 

d 2 dg 
di 

j(m — njriT 



1 + J(r'+n')m T mi+s+[2n'(me(^--^)] J 





r]2\X2-r -z—7d 2 
2n' 



2nf [ ml^ - 4- 
d 2 dg 



(154) 



The transmission function xi still has period d g /riT — 
d g /(2n'), but the first factor in the integrand has twice 
this period, dx — 2d g /n,T — d g /n'. As in the case of odd 
riT, one must choose f 2 sufficiently small to eliminate all 
but the r' = terms in the sum to ensure that the atomic 



grating has period dr- Inequalities (138) and (14C) are 
satisfied if 



2n 



which is a sufficient condition for neglect of terms with 
r' ^= in Eq. (|154|). As a result one arrives at 



f(x,t e ) 



h 



2jn'(m — n) 



s=0 



(156) 



h.(w) = f dz(l + (-lf mT 



m£-\-s-\- \a' z—w\ 



X6{[3' - [s + {a'z~ w} F )} 

a' = utO, = 2n'm£f 2 ; 

= utP = 2n j (to — n) fx; 
x 

w = — — . 



(156a) 

(156b) 
(156c) 
(156d) 



This expression can be evaluated in the same man- 



more complicated owing to two factors: (i) contributions 
h s (w) with s < [f3']j are not independent of w, and (ii) 
for s = one has to consider separately contributions 
from odd and even [a'z — w], . The situation simplifies 
for integer (3' , when the 9- factor equals 1 for s < f3' and 
for s = /?', independent of the values of z, w, and a' . 
This is the only limit considered in the subsection. For 
integer 0' , one can carry out the summation over s in Eq. 



(155) ( 156 ) 



For even f3', when X^o 1 (~ ^-T = 0' tne gratings are 
washed out and 



f(x,t) = hf 2 . 



(157) 



This result is consistent with the vanishing of the atomic 
gratings in Fig. [l(] for Talbot phases equal to n/2, 7r/4 
and 7r/6, corresponding to values of j3' equal to 2, 4, 6. 
When f}' is odd one finds 



f(x,t) = hf 2 1 + 



(-1) 



-h'(w) 



h'\ir) : : n' { 2 J dz8 



a'(3' 

1 f a' z — w 

2 " 



(158) 



(158a) 



ner used to evaluate Eq. (143), but the evaluation is where the equality 
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2 \ 2 J j 



(159) 



Eg. (|99[) with the replacements (3, a , w 
Using these values in Eqs. (101a, 101t ), one obtains 



1 CL_ W_ 

2 ' 2 > 2 



has been used. Equation ( |l58a| ) can be reduced to 



h' (w)= 2 < 



HI — 1 



{*}, 



F 



F 

w — 1 — 



< w < 2 •. 

< w < 1, 



KuK l + 2<[^} 
1 + 2 {t} f <^< 2 

< ui < 2 j^} 
2{^} f -1<«<1, 
1<w<2 {t} f . 

< w < 2 



1, 



for 



< 



for 



> 



(160) 
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FIG. 11. Gratings' amplitudes A(L) at the focal planes 
y = 2L and 7/ = |L as a function of the distance L between 
microfabricated structures. Two separated microfabricated 
structures have the same period (di = <fe) and relative width 
/i = i, /a = I- 



This expression describes an atomic grating having pe- 
riod Aw = 2 or Ax = -§■, which is n' times narrower 
than that caused by the shadow effect. For the param- 
eters chosen in Fig. [ll], values of the Talbot phase are 
equal to (fyr(m,n) = ir/n 1 (n' = 1, 3 or 5), a' = /3' = n', 
and the atomic density is given by 



f(x,t e 



ijl-i 
12 



di 



- 1 



(161) 



Since the atomic density is a periodic function of the 
distance between the MS having period Lt, the depen- 
dence of the grating amplitude A at a given focal plane 
must also be a periodic function of L having period Lt 
for a fixed value of the ratio y e /L. One period of A(L) 
is shown in Fig. [ll] at the focal planes y e = 2L, 3L/2. 
This dependence is plotted for values of L equal to ra- 
tional multiples of the Talbot length L = ^^Lt- One 
can not expect the dependence of A (L ) to be smooth, 
because the transmission function (137) is discontinuous; 
even small changes in the ratio can lead to dramatic 



changes in the atomic density (134). 



C. Comparison of the Talbot and Talbot-Lau effects 

Qualitatively, the transition from the shadow effect to 
the Talbot-Lau effect for a beam having a large angular 
divergence that is scattered by two MS parallels the tran- 
sition from spatial modulation to Talbot self-imaging of 
a collimated beam that is scattered by a single MS. Sim- 
ilarities and differences of these transitions, which oc- 
cur when the characteristic length scale in the problem 
changes from L -C Lt to L ~ Lt, can be summarized as 
follows: 
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Transition from shadow to Talbot-like profile 
(collimated beam) 



Transition from shadow effect to Talbot-Lau 
effect (divergent beam) 



Atomic density is a periodic function of the distance 
L between the MS and the screen having period Lt = 

2d 2 

Higher order gratings (with respect to the MS-grating) 
can be obtained at distances L — if, for example, 
Tlx is odd, an atomic grating having period is pro- 
duced if the MS's duty cycle / < — . 



The atomic grating's profile is the same as MS's profile, 
no compression occurs. 



Atomic density at a given focal plane y = is a 
periodic function of the distance L between the two 
separated MS, having period Lt — p(m-n) ■ 

Higher order gratings (with respect to those focused in 
the shadow-effect-regime) can be obtained at distances 
L = j^; if, for example, nj> is odd, an atomic grating 
having period -. dl is produced if the second MS's duty 
cycle / 2 < ^■ 

The atomic grating's profile is the corresponding 
shadow effect grating's profile compressed by a factor 

TIT- 



D. Additional examples including a quantum 
Talbot-Lau effect 



To make some connection with previous work, we ana- 



lyze the atomic density for the parameters of Eq. (f If ), 
corresponding to the Talbot-Lau effect studied theoreti- 
cally by Carnal et al. (1995) and realized experimentally 
by Clauser and Li (f994). The appropriate parameters 
are {m, n, j, £,],£} = {2, f, 2, f, 1,1/2}, 



d g =d u L T {2,1) 



,d 1 =2d 2 



When the distance between the MS is L 



_ L T (2,1) 



cor- 



responding to <pT = i" [the case analyzed by Carnal et 
al. (1995)], one has n T = 2n' = L T (2,1)/L = 2, and 
the corresponding values of a' and ft' obtained from Eqs. 
(156) are 



a' = 2/ 2 , = 2/l 



(163) 



As in (Carnal et al., 1995), we choose f 2 = \ and fi = \ 
or fx = ~. In order to apply the results of Sec. V.B.2, it 
is necessary that f 2 < 1/ut = 1/2; clearly, this require- 
ment is met. 

When fx = |, the par ame t er j3 ' is an integer (/?' = 1) 
and one can use Eqs. (158, 16C) to obtain the atomic 
density 



/0Me) = ~( I 



d i J F 



(164) 



coinciding with the profile obtained by Carnal et al. 
(1995). Since in this case both parameters a' and (3' 



are integers, the shadow effect does not lead to the any 
atomic grating (see below). The grating (164) arises en- 
tirely as a result of matter-wave interference. 

When fx = j , the parameter (3' — ^, and one has to 



return to Eq. ( |156| ), in which only the s = term in 
the sum contributes. Carrying out the integration in Eq. 
(156a), one finds 



0, 



for 



(162) f(x,t e ) = { 



{*} 



{*}, 



< 



for i< 



for 



< 



1 



for |< 



{ir}j 



4 
< 

< 



(165) 



< 1 



coinciding with the distribution calculated by Carnal et 
al. (1995). To compare this profile with t hat cause d 
by the shadow effect, one finds from Eqs. (102, 101c ), 
that the shadow effect distribution function is given by 
/ ^shadow = f( x ~ ^>*e). Thus, owing to matter- 
wave interference, the atomic grating (165) is shifted by 
a half-period from the grating that would have been pro- 
duced by the shadow effect alone. 

In general, the atomic gratings produced when a beam 
scatters from two, separated MS cannot be attributed 
entirely to quantum effects since the classical shadow 
effect contributes to grating formation. If the parame- 
ters are chosen in such a way, however, that the classical 
shadow effect vanishes, then any atomic gratings that 
are formed can be attributed solely to quantum matter- 
wave interference. We have already alluded to this result 



above. Returning to Eqs. (101, 102), one finds that the 
shadow effect grating S (w) disappears if the parameters 
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a = ml/2 or j3 — j(m — n)fi are integers. One can 
guarantee that a is integral by choosing 

fx = h = \ n = 1, m = 2, j = 1, I = 1, (166) 

which corresponds to y e = 2L, d± — d 2 , J ' — £ = 1, a = 
1, (3 = i and a Talbot phase (/>t = ^ Lr f 2 1) ■ The atomic 
density in this focal plane as a function of Talbot phase 
is shown in Fig. 



II 



L=7 31 Lr 



1 



L— /3|L T 






L= 7 / 31 L, 





L= I5 / 31 L r 



x/d D 



FIG. 12. An evolution of the atomic spatial distribution at 
the echo-point t — 2T (focal plane y = 2L) , induced by two 
separated microfabricated structures having the same spac- 
ings and 50% relative widths (/1 = fa = |) , with an increase 
of the distance L between structures on the scale of the Tal- 
bot distance Lt- Chosen case allows exclude any influence of 
the classical shadow effect and interpret the atomic grating 
as an exact consequence of the wave-matter interference. The 
atomic grating as a function of L is monitored for L < ^Lt 
with step ^j-I/T starting from the point where one has no in- 
fluence of the atom interference. Gr atin gs at \Lt < L < Lt 
can be obtained using the relation (13C). 



VI. TALBOT AND TALBOT-LAU EFFECTS IN A 
THERMAL ATOMIC BEAM. 

Up to this point, all effects related to a distribution 
of longitudinal velocities u in the atomic beam have 
been neglected. Averaging over u is not important for 
the shadow effect since the focal planes are located at 
y e = (m/n)L, independent of u. In both the Talbot 
and Talbot-Lau effects, the Talbot phase depends on the 
Talbot length Lt = 2<i 2 /A ( j_B, which, in turn, is pro- 
portional to u owing to the presence of the De Broglie 
wavelength. To achieve the maximum contrast in the 
Talbot and Talbot-Lau effects, it is necessary to longi- 
tudinally cool the atomic beam (Clauser and Li, 1994). 



The results of Sees. IV and VI must be averaged over 
u once changes in the Talbot phase originating from the 
distribution of longitudinal velocities becomes of order 
unity. For a thermal beam having a Maxwellian distri- 
bution over longitudinal velocities, the averaging can be 
carried out using the function tabulated by Kruse and 
Ramsey (1951). For other distributions numerical inte- 
gration is needed. Such calculations are not included in 
this contribution. 




Screen 

FIG. 13. An origin of the Talbot effect in the thermal 
beam at an asymptotic distance L > Lt- When incident 
beam of atoms having momenta p splits in a set of scattered 
beams having momenta p±nftk. The wave functions of the 
states associated with momenta p+?ik and p— ?ik, p+2?ik 
and p— 2?ik, . . . p+nRk and p— nhk acquire the same Tal- 
bot phase and therefore no dephasing between these states 
occurs, independently on the distance L and atomic veloc- 
ity u. Interference of these states remains at the asymptotic 
distance, leads to the 2nd, 4th . . . harmonics in the atomic 
distribution which form the 2— order atomic grating on the 
screen. 

Instead, we examine the role of the longitudinal ve- 
locity distribution when the width u of the longitudinal 
velocity distribution is of order of the average velocity, 



u ~ u 



(167) 



for distances 



(D/d)L T > y > L T (168) 

in the Talbot effect and separations L between the MS 

{D/d)L T > L > L T (168a) 

in the Talbot-Lau effect. We want to examine whether or 
not it is possible under these conditions to obtain atomic 
gratings having periods smaller than the MS producing 
the scattering. 
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To understand how the gratings can survive the aver- 
age over longitudinal velocities u, one should note that 
it is the atomic density ( j3(i| ) and not the wave function 
amplitude that is averaged. The phase factors in the and 
atomic density depending on the Talbot phase can be 
unity for specific combinations of the spatial harmonics 
in the atomic wave functions. The way in which this can 
be achieved is illustrated in Fig. O. When one com- 
bines on the screen components of the scattered atomic 
states associated with momenta Tik and — hk, 2hk and 
—2hk, . . . , etc., the amplitudes of the combining states 
acquire the same Talbot phases since the energy of the 
scattered atoms does not depend on the direction of scat- 
tering. The interference from these pairs of states leads 
to a superposition of atomic gratings having period d/2, 
d/4, etc.; the overall period of the grating is d/2. Grat- 
ings originated from Talbot-Lau effect can survive in a 
similar manner. 



A. Atomic density profile for a thermal beam. 



/ = £KI 2 = r^rxOO, (ma) 



d dx' 



d dx" 



x^exp[ink(2x-x' -x")]. (171b) 
The atomic density profile has a period given by 



d„ = — . 

9 2 



(172) 



Note that the density profile is independent of t for the 
times t 3> Lt/u under consideration. 



To evaluate the atomic distribution ( |171| ), it is conve- 
nient to introduce new variables 



x = - (x' + x") , x = x' - x". 



(173) 



IE.*' 



1. Talbot effect 

Consider first the Talbot effect, i. e., the atomic grat- 
ing produced when an atomic beam having negligible an- 
gular divergence, but a finite spread of longitudinal ve- 
locities, is scattered by a MS having period d. For a given 
u, the atomic density in the plane y = ut is given by [see 

Eqs. §§] f(x,t) = f- 



After adding and subtracting a term having n — in 
Eq. ( |l71b| ) one can carry out the summation to obtain 



y), making use of Eq. (|io|) . Switch- 



ing integration variables from (x',x") to (x,x) one sees 
that the 5— functions having s = and 1 are the only 
ones that contribute to the sum, and it follows from Eq. 
(|!71b[) that 



fu(x 1 1) = 2_. VriVn' ex P [* ( n -n')kx-i( 



n 2 -n' 2 ) 



; dx' 



F (x) + F [ x + - 



</>t( u )] ) 



(174) 



(169) 

where the Talbot phase 4>t{u) = Wkt = tOky/u, as given 
by Eq. (29), is a function of u for fixed y. Recall that 
r] n is a Fourier component of the amplitude transmission 
function T]{x). Since the Talbot phase is much greater 
than unity in the asymptotic region (16S), the distribu- 
tion ( |169] ) oscillates rapidly as a function of u. After av- 
eraging over it, a non-zero result arises from only those 
terms having 



n = ±n. 



(170) 



There is a constant background term / ( for a MS con- 
sisting of an array of slits, / is equal to the duty cycle / 
of the MS) corresponding to contributions with n' = n 
and an interference term f(x,t) corresponding to contri- 
butions from n' = —n. Neglecting all other terms, one 
finds 



f(x,t) = f + f(x,t), 



(171) 



where 
F(x) = 



dx I _ x \ I _ x 
—r\ \x + — \ r] \ x — — 



\x\<2 min(x,d—x) 



(175) 



For a transmission function corresponding to a peri- 
odic array of slits having duty cycle /, one finds / = 

Jo txM = /. and 



F(x) = 



2 / 




. fOT {l} F <f 



(176) 



When / < i, the two F— functions in Eq. (174) do not 
overlap with one another, and the atomic density is given 
by 

w, for < w < ^ 
/ (x,t) = f(l-f) + 2{ f- w , for | < w < f , 

0, for / < w < | 



where 



(177) 
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where w = | } . This function has period |. For 
/ > 5 one arrives at the distribution 

/ - 5, for < w < f - 5 
/(,-./)=/(] -/!-:- 2 <( ttf, " for/-i < W < I , 

f — W, for I < u; < I 

(178) 



2. Talbot-Lau effect 

To evaluate the Talbot-Lau density profile in the 
asymptotic limit (168), one must return to Eq. fll23| ) 
and average it over longitudinal velocities. Using the 
Fourier expansion of the amplitude transmission func- 



tions in Eq. (123) and setting the smooth envelope func- 
tion / (x) equal to unity, one finds 



which also has period |. The amplitude of the atomic f(x, t e ) = exp (iq [k g x — <f>r(ni, n; u) (2n 2 — mlqj^ ) 

(180) 



grating (177, 178) is given by 

A = min [/, (1 - /)] . 



,(!) 



(2) 



(179) 



X X— j(m— n)q^ n l 



(2) 



The manner in which the atomic density profile 
changes as y varies from y <C Lt to j/ > Lt is shown in 



Fig. 14 




X 

FIG. 14. Talbot effect in the thermal beam. The initial 
atomic distribution (dashed lines), created just after passing 
through the microfabricated structure with the slits' relative 
width / and spasing d, transforms into the second order gart- 
ing (having spacing |Q on the distance much larger than the 
Talbot distance owing to the wave-matter interference. 



where 4>r(rn, n;u) is given by (124). 

On averaging over u for 0y(m, n;it) 3> 1, one finds 
that only those terms in the sum having q = or 



n 2 



m£-, q^O 



(181) 



contribute to the density. The q — term is a back- 
ground, given by /1/2 [fj is defined in terms of Xj m the 
same way that / is defined in terms of \ in Eq. (171a)], 
while the terms satisfying Eq. (181) lead to the atomic 
grating. Explicitly, one finds 



f(x,t e ) = fi 



h 



dx 2 



+ f(x,t e ), (182) 



/Me) 



Q 



exp (iqk g x) x 



(i) 



F (2) 



-j(m—n)q mlq 



(182a) 



where the summation over q includes q = and other val- 



ues of q leading to integral n 2 in Eq. (y_81j). The function 
Fv is a Fourier component of the function F 2 (x) defined 
in terms of rj 2 (x) in the same way that F is defined in 
terms of 77 in Eq. (175). The density is independent of t e 
for spatial separations L 3> Lt of the MS. 

When ml is even, all q are allowed according to Eq. 
(181). When ml is odd, only even values q contribute, 
which means that the grating ( |l82a| ) has a period equal 
to dg/2. Equation ( 182a| ) has the same structure as Eq. 
(|S5"|). Repeating the calculations leading to Eq. ( |93| ) one 
arrives at the formula 



j(m— n)-l d 



f(x,te) 



j(m — n) 



E 

s=0 



dx 2 



F 2 (x 2 ) 



where 



di 



j(m — n) 



d 2 dg / J p 



1 r j even 
n , tor ml , . 

2 odd 



(183) 



(184) 



Consider now the case of MS having transmission func- 
tions (p3), when the function F 2 (x) is given by the rhs 
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of Eq. (176) with d and / replaced by d 2 and f%. Using Omitting further calculations which are essentially the 

same as those used to derive Eqs. (101, 102), one finds 



dimensionless variables 

w = £x/d gi z = x 2 // 2 d2, 



(185) 



one arrives at equations that are the analogues of Eqs. 
(p6|), namely 



/(Me) = /l/2(l-/ 2 )+/(Me), 
J2 



f(x,t e ) 



j {m-n)*- 



(186) 
(186a) 



s=0 



h s {w)= dzF(z)6((3-{s + {az-w} F )), (186b) 
Jo 



F(z) = U r Z ' f ° r f Z< ^ ! , 
' 5 - z, for z > I 



(186c) 



f2 

/■:.<•./.; - ' 2 



j (to - n) 



[(3} I +J2 F (. a s,b s ) + f (i 



8=1 



where 



F(a,b) = j dzF(z) 

■J a 



(187) 



(188) 



is given by 



2 (b 2 - a 2 ) max (a, b) < | 

F (a, 6) = ■{ 46 - 2 (a 2 + 6 2 ) - 1 a < ± < 6 

2 (6 - a) (2 - 6 - a) min (a, 6) > | 



(189) 



a = £m£f 2 , 0=j(m-n) f x . 

( F(a [a]i+1 ,l), for < w < 1 - {/3} F 

f // ^ = J F(a [a]j+1 ,l) + F(0,6 ), for 1 - {/3} F < ib < {a} F 

{W) ) F (0, b ) , for {a} F < w < 1 + {a} F - {/3} F 

[ F(G,& )+F(& Wj ,l), for 1 + {a} F - {/?} F < w < 1 



the quantities a s and fe s are given in the Appendix by 
(186d) Eqs. rtAl Ia|), and 



if {(3} F > max({a} F , 1 - {a} F ) , 

(190) 



0, 



for < w < {a} F 
for {a} F <w<l- {13} F 
F (0, b ) , for 1 - {13} F <w<l + {a} F - {(3} F 

F (0, bo) + F (b [a]i , 1) , for 1 + {a} F - {f3} F <w<l 



} , if {a} F < {13} F < 1 - {a} F , (190a) 



H > 



for < w < {a} F - {13} F 



F (a[ a ] I+ i, &h z -k, > ' i ' 

F (ffl[a] J +i, l) , for {a} F - {/3} F < w 

F (0, ho) + F (o [a]j+ i, 1) , for 1 - {/3} F < w < {a} F 

F(Q,6 ), for {a} F < ui < 1 



< 1 - {P} F 



}, if 1 - M F < {/?} F < {«} F , (190b) 



^ (fl[a] f +i)^ a ] J +i) , for < w < {a} F - {f3} F , 

i±;!f^'° lF .«{«,< !-{«),) 



0, 

F(0,& ), 



for 



for {a} F < w < 1 - {(3} F 
for 1 - {f3} F < w < 1 



(190c) 



In principle one can use Eqs. (182, 187, 19C) to derive 
a general analytical expression for the atomic density dis- 
tribution, but, given the large number of cases, such an 
expression is of limited use. Instead, one can write a 
computer code based on Eqs. (182, 187, 19C) to obtain 
the atomic density profile. Using this code, we varied 
the duty cycles of the MS to optimize the atomic grating 



amplitude in the focal planes y — ^-^-L for n = 1 — 4 and 
equal periods of the MS, j — I — 1. Calculations show 
that one has to choose 



/i — h — 2 



(191) 



in all cases except n = 3, where the optimal duty cycles 
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are given by 



Similar calculations leads to the atomic grating profiles 



f - 1 f - 1 
Ji — 2> J 2 — 4- 



(192) 



For these optimal values of the duty cycles, it is a sim- 
ple matter to obtain analytical expressions for the atomic 
density profile in a given focal plane. For example, at the 
echo plane y = 2L (n = 1, m = 2) the parameters a and 
(3 are equal 1 and i , respectively, and the atomic density 
is given by 

/Me) = l - + \{F{ ai M) 

w-^][F(0,bo) + F(b l7 l)}\, (193) 



where w 



w =F 



a i 



Using Eq. 



(189) one arrives at the atomic grating profile 



f(x,t e ) = - 



2w{l 
2w(3 



2w), for w < \ 



2w), for w > 



(194) 
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FIG. 15. Talbot-Lau effect in the thermal beam. The 
atomic density distribution at the different focal planes y 
plotted for the optimum relative widths of the separated mi- 
crofabricated structutres (/i = | and /2 = ^ in all cases ex- 
ept that for focal plane y = |L, where /2 = j). The distance 
between fields is assumed to be much larger than the Talbot 
distance. In spite on this, one can observe hi gher- orde r spa- 
tial gratings. Moreover, one finds from Eqs. ( 102 , 101) that 



for optimum widths there are no gratings caused by classi- 
cal shadow effect and conclude that using microfabricated 
structures one observes gratings originated entirely from the 
wave-matter interference. 



f ( x t ) - 1 / 9 + 2w ^ ~ 2w ^ for w< \ (194a) 

/ [x, t e ) - 36 | n _ 2w ( 3 _ 2w ), for w > \ ( ' 

at the focal plane y = |L, where w = {■jf} > 

f(x,t e )^^{ 1 + W [ 1 "l W } ^ w <\ (194b) 

J v ' eJ 8 \ 2 - w(3 - 2w), for w > \ K 1 



at the focal plane y — ^L, where w 



{*}, 



and 



f (x, te) = 



25 
27 



2w(l 
2w(3 



2w), for w < \ 
2w), for w > | 



(194c) 



at the focal plane y — |L, where w — 

atomic density profiles are shown in Fig. [l5|. 
optimal duty cycles (191, 192) correspond to the limit 
where the shadow effect vanishes, the density profiles 



. These 
Since the 



(194) cannot be vestiges of the shadow effect. They must 
originate from matter-wave interference. One can com- 
pare the density profile (194a), valid for distances L be- 
tween the MS larger than the Talbot distance, with that 
of Fig. fLOl for L ~ L T (Talbot-Lau effect). 



VII. CONCLUSION 

Atom interferometry is an emerging field of atomic, 
molecular and optical physics. In this review, we have 
focused on the scattering of atoms by one or more mi- 
crofabricated structures (MS). We have seen that the 
scattering can be described in purely classical terms for 
characteristic length scales L <C Lt, where Lt is the so- 
called Talbot length. For L > Lt, a classical description 
of the atomic center-of-mass motion is no longer ade- 
quate. Our approach has relied on an interpretation of 
the phenomena in terms of the recoil an atom acquires 
when it is scattered from a MS. With this approach, we 
could make a connection with the theory of coherent 
transients, for which a rich literature has already been 
developed. We have considered both collimated beams 
(Talbot effect) and beams having large angular diver- 
gence (shadow effect, Talbot-Lau effect) and have allowed 
for a broad distribution of longitudinal velocities in the 
atomic beam (Talbot and Talbot-Lau effects in a thermal 
beam). The next step would be to extend our consider- 
ations to regimes corresponding to Bragg scattering and 
Fraunhofer diffraction, allowing for an analysis of atom 
interferometers which split atomic wave functions into 
nonoverlapping paths. 

Scattering of atoms by MS shares both similarities and 
differences with scattering of atoms by standing- wave op- 
tical fields (SW). Similarities include a periodical recov- 
ery of the atomic interference pattern at multiples of the 
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Talbot distance (|32j) [or (|127|) for the Talbot-Lau effect], 
a compression of the atomic gratings with respect to the 
periods of the MS or SW, spatial separation of the higher 
order atomic gratings in different focal planes, and split- 
ting of the incident beam into an infinite set of scattered 
beams having momenta (p ± nhk). The differences are 
due in large part to the nature of the scattering. The 
MS produce a piecewise continuous atomic density profile 
while the SW produce a smooth atomic density profile. 
As a result, the decrease in period relative to that of the 
classical shadow effect observed in the Talbot-Lau effect 
using MS does not occur for scattering by standing-wave 
fields. Moreover, the possibility of observing a Talbot- 
Lau effect caused entirely by matter-wave interference 
(see Fig. |l2] ) does not occur for the smooth amplitude 
modulation by SW (Dubetsky and Berman, 1994). In 
the case of scattering by MS, the fact that the shadow 
effect does not give rise to atomic gratings in the focal 
plane y = 2L for MS having duty cycles fa = 4 , is di- 
rectly related to the stepwise amplitude modulation of 
the atomic beam produced by the MS, as is evident from 
Fig. 

For scattering by MS, the qualitative nature of the 
atomic density profile depends on the properties of the 
incident atomic beam. When one observes the Talbot 
effect using a monovelocity beam, the atomic gratings 
are discontinuous functions (see Fig. |||), but when one 
averages these gratings over longitudinal velocities, the 
atomic density is transformed into a piecewise continu- 



ous profile (compare with Fig. 14). Similarly, the shadow 
effect and Talbot-Lau effect atomic gratings (which are 
averaged over transverse velocities) are piecewise contin- 
uous, but they are transformed into profiles which are 
discontinuous only in the second derivative when aver- 
aged over longitudinal velocities (compare Figs. [l(] and 
|l5| ) . These examples show that averaging over transverse 
or longitudinal velocities tends to smooth out the atomic 
density profiles. 

It is clear that many of the situations analyzed in this 
chapter have direct applications to atom lithography. We 
can expect that future developments in this emerging 
field will incorporate many of the basic ideas which have 
been encountered in our discussion. 
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APPENDIX: 



A method for calculating the integral 



h m (w) = / dz9 [{/3} F - {az - w} F ] , (Al) 
Jo 

with < w < 1 is presented here. The function 9 is the 
Heaviside step function and [a]i and {a}i? refer to the 
integral and fractional parts of a. Depending on the val- 
ues of the parameters a, (3 and w the integrand in this 
equation can jump from to 1 several times inside the 
range z £ [0, 1] . The value of the integral is the total 
length of the intervals in this range where 



C (z) = {az - w} F < {f3} } 



(A2) 



To determine this length one needs to find the values of 
z for which the function {az — w} F equals and where 
it equals {[3} F .This function is shown in Fig. 16. 



{ocz-w} F 



z=l 




b aj bj 



a [al b [aj a M,-H b M , 



FIG. 16. Function {az — w} undergoes jumps at the 
points a r — (j — 1 + w) /a and is equal to {13} at the points 
b r = (r - 1 + w + {13}) /a. For H > 1 and 1 < r < [ a ], 
a r < 1, i. e. contribution to Eq. from interim [ai,a[ a j] 
is ([a] — 1) (b r — Or); when w > 1 — {(3} one gets term bo from 
[0, ai]; term from z £ [a[ a ],6[ a ]] is equal to {(3} /a or l-a^] 
for b[ a ] < 1 or 6[ a i > 1; in the same manner one calculates 
term from interim [a[ a j +1 , 6[ a ] + i] in different possible cases. 



When az — w = r — 1, where r is an integer, £ (z) = 0, 
i. e. zeros of £ (z) are given by 



z = a r = (r — 1 + w) /a. 
For z > a r the function ( (z) evolves as 



(A3) 
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( (z) — a (z — a r ) 
and equals {f3} F at the point 
{(3} F 



(A4) 



a r 



(r-l + w + {P} F )/a. (A5) 



Since contributions to the integral (|Al| ) vanish unless 
a r < 1 and b r > and since < w < 1, it follows 
from Eqs. (|a|, |a|) that < r < (1 + [a]/). For the time 
being, we assume that a > 1. 

All intervals [a r ,b r ] totally or partially within the 
range [0,1] contribute to the integral (Al). Let us de- 
note the contribution from the range [a r , b r ] by A r and 
the total value of the integral by 



i+H 



VhH = Ar - 



(A6) 



r=0 



When a r > and b r < 1 the interval [a r , b r ] lies entirely 
to the range [0, 1] and 



. 1 , 



a r > and b r < 1 (A7) 



When 1 > a r > and b r > 1, the maximum value of z 
contributing to the integral ( |A1[ ) is z = 1 and 

A r = 1 - a r \ 1 > a r > and o r > 1 (A8) 

Similarly, for a r < and < o r < 1, 

A r = 6 r ; a r < and < b r < 1 (A8a) 

and for 6 r < 0, 

A r = 0; 6 r < (A8b) 

For given r, a, and /?, the values of a r and b r can depend 
on to, giving rise to a dependence of htg], on w. Note 



that (b r — a r ) < 1, which follows from Eqs. ( A5) and the 
assumption that a > 1. 

We first consider the range 1 < r < [a] 1 — 1 , for which 
b r < 1 and a r > for any w £ [0, 1]. It then follows from 
Eq. (A7) that the total contribution A' to the integral 
(Al) from the region [ai, 6[ a ] _i] is given by 



A' = J2 A r ^ {P)f (Hz " !)/«■ 



(A9) 



This contribution is independent of to and represents a 
constant background term. 

Since r < [a]j + 1, the only remaining contributions 
to the integral can come from Aq, and A[ a ] +i. 

These terms depend on w and represent the atomic 
gratings. Let us first consider Aq. If r = 0, bo = 
(-1 + to + {/3} F )/a < 1 and a = (- 1 + w)/a < 0. 
It then follows from Eqs. |A84 |A8b| ) that 



(to + {[3} F ~ I) /a, for to > (1 - {/3} F ) = to 







for to < 



(A10) 



We now turn our attention to j4[ q j i and ^4[ Q ] +i. It 
follows from Eq. (A3) that dUj G [0, 1]. Only the points 
\a\ji a [a] 7 +ii b[a] +i can lie to the right of the range 
[0, 1], which occurs when 

to> ({a} F + l- {f3} F ) = w l7 

w>{a} F = w 2 , (All) 

w > ({«}f - {/9} F ) = w 3 , 



respectively. From Eqs. ( |A8| , A8a, A8b), one finds that 
A[ a ]j, and A [a]i+1 are given by 



A, 



_ f b Hi - a [a], = {/3}f/" w < tOl 



1 — Ol 



a]i 



{a} F —w + l 



W > tOl 



(A12) 



.4 



Hj + 1 



to > W2 

1 - a [a]j + l = W 3 <W <W 2 

b[ a ]t+i - a[ a ] 1+ i = {/3} F /a w < to 3 

(A12a) 

In order to sum Ao, A[ a ^, and j4[„] +1 it is convenient 
to separate regions of a and (3 according to the relative 
values of toi, t02, w 3 and to. Since 



W3 < w 2 < tOi, 



(A13) 



u>3 < to < toi, (Al3a) 

and toa = {q;}^ > one can distinguish four cases 

w 3 < < w < w 2 ; [{/3} F > max({a} F , 1 - {a} F )] , 

(A14) 

w 3 < < w 2 < w; [{a} F < {/3} F < 1 - {a} F } , 

(Al4a) 

0<w 3 <w<w 2 ; [1 - {a} F < {I3} F < {a} F } , 

(A14b) 

< w 3 < w 2 < w ; [{(3} F < min ({a} F , 1 - {a} F )} . 

(Al4c) 



Consider, for example, the case ( A14 ) for the range of to 
given by 



< to < to. 



(A15) 
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For integer (3 this corresponds to the entire range of al- 
lowed w, < w < 1. From Eqs. jAlOj |A12| , |Al2a 
one finds A = 0, A [a]j = {{3} F /a, and A [a]i+1 = 
({a} F — w) I a. As a result, one finds that A (w) = 
+ A [a]i + A [a]j+1 is given by 



A(w) = {{/3} F + {a} F - w) /a. 



By combining Eqs. ( |A9| , A16), one obtains 



(A16) 



(A17) 



Even though Eq. ( |A17| ) have been derived for a > 1, one 
can verify that it holds for arbitrary a. 



Other values of w and other cases ( A14a - A14c ) can be 
considered in the same manner. As a result, one arrives 



at Eqs. (|100|, |l0lj) of the text. 
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[1] The scattering of atoms by standing-wave fields, rather 
than MS, is a bit more subtle. For resonant standing wave 
fields, which can act as amplitude gratings, the situation 
is unchanged. On the other hand, off-resonant fields act 
as phase gratings for the atoms; as such, they produce no 
effect on classically moving particles. Strictly speaking, 
therefore, one must quantize the center-of-mass motion to 
calculate the scattering of the atoms by the fields. Never- 
theless, in a manner analogous to the normal photon echo, 
it is possible to assign phases to the atoms while they are 
freely evolving between the MS and the screen and still 
consider the motion as classical in these regions. 
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